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Dear  Sir/Madam: 

This  letter  constitutes  the  revised  final  report  for  ONR  contract  N00014-88-K-0203 
awarded  to  Old  Dominion  University  for  the  period  October  1,  1987  though  January 
31,  1989.  The  original  final  report  was  submitted  to  the  Office  of  Naval  Research  as 
part  of  a  renewal  proposal. 

The  bulk  of  the  research  activities  conducted  under  this  contract  consisted  of  basic 
research  on  ocean  flow  dynamics  as  it  pertains  to  the  prediction  of  ocean  motion. 
During  the  contract  period,  three  papers  describing  aspects  of  the  research  were 
published  in  the  adjudicated  scientific  literature.  These  papers  are: 

1.  “Genesis  of  the  Gulf  of  Mexico  Ring  as  Determined  from  Kinematic  Analysis,” 
J.  Geophys.  Res.,  92(C11),  11727-11740,  1987. 

2.  “Observed  and  Simulated  Kinematic  Properties  of  Loop  Current  Rings,”  J.  Geo¬ 
phys.  Res.,  93(C2),  1189-1198,  1988. 

3.  “Notes  on  the  Cluster  Method  for  Interpreting  Relative  Motions,”  J.  Geo¬ 
phys.  Res.,  93(C8),  9337-9339,  1988. 

In  addition  to  these  papers,  the  results  of  the  research  were  reported  at  the  1988  Liege 
Colloquium  on  Ocean  Hydrodynamics,  Mesoscale/Synoptic  Coherence  in  Geophysical 
Turbulence.. 

As  a  part  of  the  research  effort,  the  funds  were  used  to  support  William  Iiidest,  a 
Ph.D.  graduate  student  in  the  Department  of  Oceanography. 

I  am  grateful  to  the  Office  of  Naval  Research  for  continued  support  for  this  research. 
Respectfully, 

A.  D.  Kirwan,  Jr. 

Samuel  L.  and  Fay  M.  Slover 
Chair  of  Oceanography 
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Genesis  of# a  Gulf  of  Mexico  Ring  as  Determined  From  Kinematic  Analyses 

James  K.  Lewis 


ScifiKv  ApplUuliiim  htlfnialiniml  Corptiriilion,  CoUeqe  SlalUm,  Texas 

A.  D.  Kirwan,  Jr.’  < 


Umrersiiy  iif  Simili  Tturula.  Oeparinu'M  of  Marine  Seiciice.  St.  Petersburg 

The  kineniiiiics  of  the  Loop  Current  are  siudieti  using  trajectories  of  drirters  in  the  Gulf  of  Mexico 
during  mid-June  through  September  1985.  One  of  the  drifters  was  in  the  Loop  Current  proper,  while 
tilhur  (Iriftvrs  were  in  twii  recently  shed  l.iHip  ('urrunt  rings.  The  drifter  in  the  Loop  Current  showed 
strong  aniieyclonic  motion  during  the  study  period.  This  Loup  Current  anticyclone  Hrst  began  off  the 
northwest  coast  of  Cuba.  It  rapidly  moved  northward  into  the  Gulf  of  Mexico  as  a  ring  pinched  off  from 
the  Loop  Current.  Analysis  of  the  Loop  Current  drifter  motion  showed  that  the  anticyclone  became  an 
integr.il  part  of  the  Loop  Current,  taking  on  many  of  the  characteristics  of  the  most  recently  shed  ring. 
The  results  of  the  .inalysts  suggest  .1  process  by  which  Loop  Current  nngs  can  be  generated.  Apparently, 
tins  mecliatiism  cati  cause  the  Loop  Current  to  become  reconligured  in  2-3  months  for  beginning  the 
process  of  ring  separation. 


I.  iNTRODt'CTtON 

The  shedding  of  Loop  Current  ring.s  has  a  major  impact  on 
processes  in  the  central  and  western  Gulf  of  Mexico.  Thc.se 
liirgc  anticyclones  transport  a  tremendous  amount  of  momen¬ 
tum.  heat,  and  salt  across  the  gulf,  all  the  way  to  the  Mexican 
coiist  ir.lliot.  I982:  Kinmi  et  al.,  I984;  Lewis  and  Ktnsan. 
I98SJ,  In  Older  to  consider  hahinccs  of  inomcntiim.  mass,  and 
heat  within  the  Gulf  of  Mexico,  it  is  important  to  have  some 
idea  of  the  charactciisiics  of  the  kinematics  of  the  Loop  Cur¬ 
rent.  including  how  often  a  ring  may  be  pinched  off. 

During  June  1 985.  an  attempt  was  made  lo  put  an  Argos 
drifter  into  the  Loop  Current  as  a  ring  was  pinching  off.  How¬ 
ever.  the  ring  did  not  totally  discoiincct  from  the  Loop  Cur¬ 
rent.  and  the  driller  ended  up  spending  approximately  3 
months  in  the  Loop  Current  proper.  The  drifter  exited  the 
Gulf  of  Mexico  (GOM)  through  the  Florida  Straits  in  Septem¬ 
ber  1985  (Figure  I,  drifter  .7.354).  Fortunately,  a  second  drifter 
(.3378.  Figure  I)  was  placed  in  the  ring  in  July  1985.  Together, 
these  two  Lagrangian  data  sets  provide  a  unique  view  of  the 
motion  of  the  Loop  Current  while  it  is  tn  the  process  of  shed¬ 
ding  a  ring.  They  provide  special  insight  into  Loop  Current 
kinemalics  as  this  system  extends  northward,  as  well  as  the 
kinematics  of  a  Loop  Current  ring.  Finally,  wc  comment  on 
the  existence  of  a  third  anticyclone  in  the  GOM.  The  eddv 
was  di.scovered  serendipitou.sly  as  a  third  drifter  became  en¬ 
trained  in  its  (low  field  and  circled  around  the  eddy  all  the 
w.ty  to  the  Mexicar,  coast  (G.  Forristall.  personal  communi¬ 
cation.  1986),  This  eddy  hits  been  referretl  to  as  Ghost  Eddy 
because  there  had  been  no  previous  indication  of  its  existence 
prior  lo  the  first  of  August  1985.  Its  general  motion  and  extent 
during  August  and  September  1985  is  shown  in  Figure  1. 
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The  following  analysis  is  obtained  from  the  movement  of 
drifters  3354  and  3378  as  well  as  concurrent  sea  surface  tem¬ 
perature  (SST)  data  and  XBT  data.  These,  along  with  the 
presence  and  location  of  Ghost  Eddy,  suggest  a  new  process 
by  which  an  anticyclonic  vortex  is  formed  that  can  eventually 
become  a  Loop  Current  ring.  The  original  mechanism  for 
spinning  up  the  vortex  appears  to  be  the  lateral  sharing 
stresses  cau.sed  by  the  Loop  Current  off  the  northwest  coast  of 
Cuoji. 

A  detailed  description  of  the  kinematic  characteristics  of  the 
l.oop  Current  and  the  ring  is  provided  by  an  analysis  of  the 
ir.uer  paths.  This  analysis  gives  us  the  time  histories  of  the 
rotation  rate,  eccentricity  and  orientation  of  the  ellipses  of  the 
trajectories,  swirl  velocities,  and  movement  of  the  centers  of 
rotation.  A  comparison  is  made  of  the  kinematics  of  the  Loop 
Current  as  ceterniined  by  drifter  3354  during  the  time  it  co¬ 
existed  with  a  ring.  The  latter  is  designated  by  drifter  3378. 
The  comparison  of  the  last  part  of  drifter  3354  (before  it  exited 
the  GOM)  with  the  first  month  of  ring  3378  provides  an  in¬ 
teresting  contrast  between  the  motion  characteristics  of  a 
northern  extension  of  the  Loop  Current  before  and  after  it 
pinches  off  to  become  a  ring. 

2.  Path  Data 

The  path  data  used  in  this  study  are  the  position  data  of  the 
drifters  with  Argos  identification  numbers  3354  and  3378.  The 
GOM  ring  that  was  seeded  by  drifter  3378  will  be  referred  to 
as  ring  3378.  References  to  drifter  3354  will  indicate  a  refer¬ 
ence  to  the  Loop  Current  proper.  Drifter  3378  was  drogued  by 
a  weighted  200-m  line,  while  drifter  3354  was  drogued  by  a 
window  shade  drogue  at  l(X)m. 

Drifter  3354  was  seeded  at  25.9'’N.  87.9'’W  on  June  18,  1985. 
The  drifter  immediately  moved  southeastward  ~525  km.  and 
reached  2.3.2''N.  83.7"W  by  June  29,  1985.  At  this  point  the 
drifter  became  entrained  in  a  westward,  anticyclonic  flow  field 
with  a  center  of  rotation  at  about  24°N,  85.5°W,  After  two 
rotations  (into  mid-August)  (he  drifter  suddenly  moved  north¬ 
westward  and  made  three  additional  rotations  centered  at 
about  25.5' N.  86.5''W.  The  drifter  then  left  this  flow  pattern  in 
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Fig.  I.  Traj«clorie<  of  drifler  .V'54  from  mitl-June  Ihrough  mid-Sepiemh«r  I98S  and  of  drifler  3378  from  mid-July 
through  Seplemhcr  198.5.  Also  shown  is  the  locnliun  of  Ghost  Eddy  duiing  August  and  September  198S.  Squares  denote 
the  beginning  povilnins  of  the  onflcr  ir.njeclories  and  triangles  denote  the  end  positions. 


niid-Scpiember  1985  and  exiled  the  GOM  through  the  Flori¬ 
da  Straits. 

Drifter  3.378  was  seeded  in  ring  3378  at  26.4  N.  89.3  W  on 
July  18.  1985  IJy  that  time,  the  ring  had  totally  scptirated 
from  the  Loop  Giirreni  As  shown  in  l•■ipure  I.  ring  3.378 
slowly  moved  weslw.ird.  retichiiig  approxiniiitely  the  91  W 
meridian  by  mid-August  Recall  thiit  after  mid-August,  the 
Loop  Current  (as  determined  by  drifter  3354)  moved  north¬ 


westward  and  remained  there  until  at  least  mid-September. 
During  that  .same  30-day  period,  ring  3378  continued  moving 
westward,  reaching  to  the  92.5' W  meridian. 

The  ollshorc  industry  launched  a  drifter  in  June  1985  which 
eventually  became  entrained  in  the  flow  field  of  Ghost  Eddy. 
At  the  beginning  of  August  1985.  this  ring  was  centered  at 
2.3.5  N.  93  W  (I'igure  I).  The  ring  gradually  moved  southwe.st- 
ward  across  the  deepest  portion  of  the  GOM.  By  the  end  of 
September  1985.  Ghost  Eddy  was  still  strongly  rotating  at 
23  N.  94.5  W. 


FABI.r  t  l.isling  of  Criiives  During  Winch  .XHT  Dala  Were 
Cullecled  in  the  Ccnir.il  and  Easiern  Gulf  «f  Mexico 


Vessel 

Diilcll9S.S) 

F..  M.  Qiimiv 

May  7  8 

M.  V  Ni-suir  1 

Mav  It)  II 

M  V  Siena  llispaimi 

May  17  18 

K.  M.  Qiiivnv 

Mav  26  27 

M  'V  Siena  Hiipaina 

May  21  28 

M'V  Siena  lli^pania 

May  .ttJ-31 

E  .\f.  Qneenv 

June  7  8 

L.  M.  Qaeenv 

June  13-14 

M.'V  Siena  lliipania 

June  26-28 

M/V  Siena  lli.fpania 

June  29-.30 

t.  Af.  Qiieeny 

July  1-2 

R  'V  Siincaasicr 

July  9  22 

M  V  Nal  Co  6 

July  16-19 

M'V  Siena  Hhpuma 

July  16-17 

E.  M.  Qiieenv 

July  31  to  Aug.  1 

M'V  Siena  lliipania 

Aug.  16-18 

M  V  Anihaiiailor 

Sepl.  4  3 

M.  V  Anihaiiadiir 

Sepi,  13  14 

M'V  iiiih/ittiithir 

Scni  24-2.5 

3.  ThMFhRATURK  DaTA 

A  number  of  XBT  dala  sets  were  collected  in  the  easiern 
GOM  during  May-Seplembcr  1985  (Table  I).  These  data  pro¬ 
vide  indications  of  the  location  and  the  structure  of  the  Loop 
Current  and  ring  3378.  In  addition,  weekly  SST  contour 
charts  arc  available  for  the  entire  gulf.  Some  of  these  SST  data 
arc  prc.scnlcd  here,  but  unfortunately  little  in  the  way  of  ther¬ 
mal  structure  can  be  picked  out  .since  the  SST  gradients  are 
small  m  the  GOM  during  the  study  period  (summer  and  early 
fall). 

Cruise  tracks  for  the  period  May  26-31.  1985.  arc  shown  in 
Figure  2a.  The  resulting  XBT  data  are  shown  in  Figures 
2/»  -2(1.  These  datti  were  collected  before  ring  3378  was  shed, 
and  the  edges  of  the  Loop  Current  (maximum  horizontal  tem¬ 
perature  gradients)  indicated  by  the  XBT  data  are  shown  in 
Figure  2a.  (Smaller  slopes  of  the  isotherms  were  taken  as  an 
indication  that  the  transects  were  crossing  the  Loop  Current 
edge  at  an  angle,  and  some  of  the  edges  were  thus  drawn  at  a 
45  .'ii'.gle  to  the  cruise  track.)  it  appears  that  the  Loop  Cur- 
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Fig.  2.  (a)  Cruise  irseks  for  XBT  data  culiccied  during  May  26-31.  1985.  The  arrows  denote  the  flow  at  the  edges  of 
the  Loop  Current  based  on  the  vertical  temperature  structure  shown  in  (h)  temperature  data  from  the  £.  M.  Quteny  cruise. 
May  26-27. 198.4.  (,.)  temperature  data  from  the  Stenn  >lt.%pamii  cruise.  May  27-28.  1985.  and  (d)  temperature  data  from  the 
Siena  Uispama  crui.se.  May  3(1  31. 1985. 


rent  was  flowing  northward  on  the  east  side  of  the  Yucatan 
Straits,  swirled  westward  and  around  the  region  where  ring 
.3378  would  be  spawned,  and  then  turned  northeastward 
before  exiting  the  COM. 

A.s  previously  mentioned,  drifter  3354  was  deployed  in  the 
northward  extension  of  the  Loop  Current  during  mid-June 
1985.  The  trajectory  of  3354  during  June  18-24.  1985.  is  plot¬ 
ted  on  the  corresponding  COM  SST  chart  in  Figure  3.  As 


indicated  by  the  trajectory  (and  to  some  degree  by  the  SST 
data),  ring  3378  had  obviously  not  separated  from  the  Loop 
Current  during  this  period.  The  return  flow  from  the  area  of 
the  northward  extension  appears  to  be  mostly  southeastward, 
toward  the  northern  coast  of  Cuba. 

Three  XBT  data  sets  were  collected  during  June  26  to  July 
2,  1985  (Figure  4),  and  the  temperature  profiles  are  shown  in 
Figures  4h-4d.  The  outline  of  ring  3378  is  quite  distinct  in 
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Tig  3.  1  filjeclory  of  drifter  3354  a;id  SST  data  for  June  18-24. 1985. 


I'lgtircs  4/i  iind  4c.  while  liic  Loop  Current  can  be  seen  to 
extend  to  25  N  in  Figure  4il.  These  data  would  seem  to  imply 
tli.it  ring  .3378  is  still  connected  to  the  Loop  Current.  The  SST 
data  fv'r  June  29  to  July  5.  1985.  along  with  the  associated 
trajectory  of  drifter  3354  arc  shown  in  Figure  5.  The  outlines 


of  the  ring  and  the  Loop  Current  from  the  XBT  data  in 
Figure  4ii  arc  also  shown,  and  -ve  see  some  of  the  rotational 
characteristics  within  the  Loop  Current  from  the  path  of  drif¬ 
ter  3354.  The  rotation  near  the  Cuban  coast  continued  into 
mid-July. 


Fig.  4.  (d)  Crtiisu  tracks  fur  Xlif  data  collected  during  June  26  to  July  2.  1985.  Hie  arrows  denote  the  flow  at  the 
edges  Ilf  ilie  l.uiip  Ctirreiii  based  on  the  vertical  icnipcraliirc  structure  shown  in  (h)  temperature  data  from  the  Steiui 
//ifpiimii  cruise.  Jui.e  26-28.  1985.  (r|  temperature  data  from  the  Steiia  Ulsponia  cruise,  June  29-30,  1983,  and  (<0 
Iciiipcralurc  data  from  the  I..  M.  (Jinriiv  cruise.  July  1-2.  198.5. 
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During  July  16  19.  1985.  two  XBT  diila  sets  were  collected 
iFigure  6)  wliich  distinctly  show  that  ring  .3.378  luid  separated 
from  the  Loop  Current.  The  .SST  contour.s  for  the  same  period 
are  shown  in  Figure  7.  There  is  no  surface  signature  of  ring 
.3.378,  although  its  rotation  is  well  delineated  by  the  trajectory 
of  drifter  .3378.  As  for  the  Loop  Current.  Figure  7  shows  that 
the  rotational  feature  within  the  current  had  reached  the 
northern  edge  of  the  current,  about  25'  N. 

Both  the  Loop  Current  aiul  ring  .3.378  coiitiiuicd  to  rotate 
anticNclonically  through  August  1985.  An  XBT  data  set  col¬ 
lected  on  August  16-18.  1985  (Figure  8)  indicates  that  the 
Loop  Current  had  e.xtendcd  northward  to  26.6  N,  Looking  at 
the  SST  map  for  August  12  -19,  1985  (Figure  9).  we  .sec  that 
ring  3.378  and  the  Loop  Current  were  rotating  at  about  the 
same  hitiiiide.  25.5  N  By  nud-Scpieiiilier.  the  sea  surface  had 
cooled  sufficiently  in  the  eastern  (iOM  so  that  the  northern 


extension  of  the  Loop  Current  is  readily  defined  (Figure  10). 
The  trajectories  for  this  period  of  time  show  ring  3378 
strongly  rotating  but  drifter  3354  leaving  the  ro'ational  fea¬ 
ture  of  the  Loop  Current  along  the  eastern  side  of  a  28‘’C 
tongue  of  water.  The  following  week's  SST  chart  (Figure  11) 
shows  the  28  C  tongue  reaching  up  to  27.5°N,  while  driller 
33.^4  moved  toward  the  Florida  Straits  along  the  northern 
coast  of  Cuba.  Ring  3378  was  still  rotating  with  a  center  at 
about  25.5  N.  92  W. 

4.  KisKsiATic  Analysis 

The  trajectory  data  were  used  to  calculate  various  kin¬ 
ematic  paramcicrs  of  the  (low  field.  The  methodology  outlined 
by  Kinwn  et  al.  [1984,  I98T)  for  geophysical  fluids  was  used 
in  these  calculations.  (I’or  more  details,  see  Kirwan  et  al. 
{ 1987]  )  Although  Kirwdii  et  al.  [1984]  used  a  l(X)-hour  low- 
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Fig.  5.  Trajectory  of  drifter  3354  (large  arrowi  and  SST  data  (  C)  for  June  29  to  July  5,  1915.  Shorter  arrows  denote 
the  flow  at  the  edges  of  the  Loop  Current  while  dolled  lines  indicsle  the  locations  of  Loop  Cuntsil  waters  (based  on  XBT 
data) 
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Tig.  6.  (ill  Crune  tracks  for  XBT  data  collected  during  July  16-19.  1985.  The  arrows  denote  the  flow  at  the  edges  of 
the  l.uop  Current  and  ring  based  on  the  vertical  temperature  structure  shown  in  (ft)  and  (r|  both  temperature  data  from 
the  M,  V  Nat  Co  6  cruise.  July  16-19.  1985.  and  (d)  temperature  data  from  iheSli'itu  JJispnma  cruise,  July  16-17,  1985. 


pass  liltcr  in  making  their  calculation '•  for  a  ring  in  the  west¬ 
ern  COM.  It  was  found  in  this  study  that  such  a  liltcr  was  not 
sufficient  to  remove  higher-frequency  fluctuations  recorded  by 
drifter  .J.7.54.  As  a  result,  we  used  a  164-hour  (half  power  point) 
low -pass  filler  on  the  velocity  data  of  the  drifters.  We  then 
calculated  the  period  of  rotation,  the  ellipticity,  the  orienta¬ 
tion.  and  the  velocity  about  the  translation  ring  center  (swirl 
velocity)  as  seen  by  drifters  .J354  and  3378. 

hmp  Current  Kiiwittalics 

The  ob.served  (fillcrcd)  Loop  Current  velocities  as  seen  by 
drifter  3354  arc  shown  in  Figure  12.  The  speed  fluctuates  in 
magnitude  from  15  to  KO  cm/s.  and  there  are  minimums  at  the 


beginning  of  tiie  record  and  during  Julian  days  210-225.  It 
was  during  days  220-230  that  drifter  3354  moved  northwest¬ 
ward.  an  apparent  rapid  northward  extension  of  the  Loop 
Current.  The  swirl  velocity  associated  with  the  rotation  of  the 
Loop  Current  is  shown  in  Figure  13.  These  magnitudes  also 
vary  from  1 5  to  80  cm/s. 

The  extent  of  water  involved  with  the  northward  movement 
of  the  Loop  Current  is  indicated  by  Figure  14,  which  shows 
the  time  history  of  the  distances  from  the  drifter  to  the  center 
of  rotation.  Prior  to  the  northwest  movement,  the  maximum 
radius  of  the  rotation  as  seen  by  the  drifter  was  85  km.  After 
the  northwest  movement  the  radius  increased  up  to  100  km. 

The  rotational  frequency  of  the  Loop  Current  is  shown  in 
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Tig  7  Tr.'ijecloricA  of  drifters  3354  (solid  Ijrye  arrow  I  and  3378  (dashed  large  arrttw)  and  SST  data  (  C)  for  July  16-23, 
l‘(X5.  Shorter  arrows  denote  the  dow'  at  the  edges  of  the  Limp  Current,  while  dotted  lines  indicate  the  locations  of  Loop 
Current  waters thased  on  XBT  data) 

higtire  15.  (hor  the  calculation  of  this  term,  .sec  Kiruum  el  al.  mosed  northwest,  the  mean  period  of  rotation  decreased  from 
ri9K4]  I  Figure  15  shows  it  general  decrease  and  then  increase  about  1.3  days  to  approximately  105  days.  If  we  compare 
in  frequency,  with  a  minimum  occurring  right  before  the  Figures  14  and  15.  we  see  that  longer  periods  tend  to  be 
iiotthwtr-twatd  cxicnston  of  the  Witter  mass.  As  the  water  asstKiated  with  larger  radii. 


fig  8  to)  Cruise  tr.ick  for  XBT  data  collected  during  August  16- 18.  1985.  The  arrow  denotes  the  flow  at  the  edge  of 
the  Loop  Current  hiKcd  on  the  vertical  tcinperatiire  structure  shown  tn  Ih)  temperature  data  from  the  Slemt  Hispanta 
cruise.  August  16  18.  198.4 


I1.7.U 


t  rwis  AND  KmwAS  Ciniximit  a  Ciuir  <ir  Mrxn  o  Rinci 


I'ig.  9.  TrajccloricA  (if  drifters  (Ai'lid  large  arrow)  and  JJ7R  (daahed  large  arrow)  and  SST  data  ('C|  for  August 
12  19.  I9R.^  .Shiirier  arrow  denoles  llie  How  at  the  edge  of  the  Loop  Current,  while  doited  lines  indicate  the  location  of 
Loi.'p  C'urreiil  walcrb  (hased  on  XDT  dulal. 


As  Tor  the  shape  of  the  lield  of  rotation,  we  use  eccentricity 
r.  defined  as  the  major  axis  leiigih  divided  by  the  minor  axis 
length  Also,  we  calciihiic  the  north. ctisl  orientation  of  the 
m.ijor  axis,  and  both  of  these  v.iritiblcs  arc  shown  in  I'igure 
16  The  How  lield  starts  out  rather  elliptical  l»*  -  1.7.^)  but 
becomes  more  circtil.ir  b\  da\  25(1  (e  ■=  I  4l.  The  elliptical 
orientation  was  mostly  cast  west  at  lirst  but  eventually 
became  more  northwest  souihcast. 


Kiiicimtlics  of  Riiui  JS7S 

The  filtered  velocities  of  ring  3.’78  are  shown  in  Figure  '7. 
The  magnitudes  of  tliese  oscillations  are  relatively  large,  up  to 
~85  cin.'s  with  a  minimum  of  about  SO  cm/s.  The  swirl  speeds 
are  shown  in  Figure  18.  and  we  see  here  an  initial  increase 
from  20  to  75  cm/s  followed  by  a  decrease  to  about  50  to  60 
cm/s.  These  variations  in  swirl  magnitude  coindde  with  vari- 


hia  It)  Trajaloncs  of  drifters  .t.1.S4  |solid  large  arrow)  and  .U78  (dashed  large  arrow)  and  SST  data  ('Q  for  September 

10  17. 1985. 
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Fig.  1 1.  Tnijeciories  of  drifters  .1354  (solid  large  arrow)  and  3378  (dashed  large  arrow)  and  SST  data  ('  C)  for  September 

17-24.  19*5. 


attons  in  the  size  of  the  circle  of  rotation  (Figure  I).  Figure  19 
quantities  the  distances  from  the  center  of  rotation,  with  a 
nia.ximutn  radius  of  95  km. 

The  rotation  frequency  of  ring  3378  is  shown  in  Figure  .’O. 
1  he  initial  minimum  in  frequency  (~0.l  1  ')  is  followed  by  .m 
increase  to  about  0.14'*.  From  about  day  220  and  on.  the 
rotational  frequency  of  ring  3378  gradually  dccrca.sed.  with  the 
period  of  rotation  increasing  from  about  8  days  to  about  9 
tlays. 


Loop  Currem  (33541 


Julian  Days  (1985) 

Fig.  13.  Time  histories  of  the  filtered  observed  speed  comiwnenls 
for  drifter  .'354.  The  light  curve  is  the  east/wesl  speed,  and  the  darker 
curve  IS  the  north, 'soulh  speed. 


Considering  the  shape,  the  data  indicate  that  the  ring  was 
initially  rather  elliptical  (Figure  21).  Within  two  rotation 
periods  this  eccentricity  decreased  somewhat,  with  the  ring 
becoming  almost  circular.  The  results  shown  in  Figure  21 
imply  that  the  inilipl  ellipse  orientation  was  mostly  cast/west. 
This  orientation  did  not  seem  to  change  significantly  as  the 
ring  became  more  circular. 

5.  Discussion 

flic  cliissiciil  concept  of  Loop  Current  processes  is  one  in 
which  part  of  the  Gulf  Stream  at  first  flows  directly  from  the 
Yucatan  Straits  to  the  Florida  Straits.  Within  this  flow  field. 


|  jp.  13  Time  hisiories  of  the  swirl  speed  components  about  the 
center  of  rpiation  for  the  Loop  Current.  The  light  curve  it  the  eait/ 
wevi  s|iccd,  and  (he  darker  curve  it  ihc  norih/soulh  ipeed. 
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Kig.  14  Time  hislories  of  Ihe  east/w;st  (light  curve)  and  north/ 
Aouth  tdarker  curve)  divtatices  of  drifter  3. 54  from  the  center  of  rota¬ 
tion  in  the  Loop  Curreni 


Julian  Days  (19BS) 

Fig.  16.  Time  hislories  of  the  Loop  Curreni  eccentricity 
(asterisks)  and  orientation  of  major  axis  (dashed  tines).  For  the  flow 
field  orientation,  north  is  in  the  positive  t  directio.i,  and  east  is  in  the 
positive  lime  direction. 


inslahilitic.'i  exist  which  result  in  meandering  of  the  Loop  Cur¬ 
rent  lllwibiiri  ami  Thompson,  1980].  As  the  states  of  the  mean¬ 
ders  increase  and  reach  northward,  it  is  generally  believed  that 
the  flow  (ielvl  wraps  back  onto  itself  and  “shorts'*  across  the 
stream  of  (low  •  part  of  the  How  would  still  go  northward 
around  the  Loop  Current  extension,  while  the  remainder  of 
the  How  would  take  the  more  direct,  .souliierly  route  to  the 
Florida  Straits  This  is  .inalogous  to  other  geophysical  phe¬ 
nomena.  such  as  the  creation  of  ox-bow  lakes  by  meandering 
rivers,  (■inally.  as  the  curvature  of  the  (low  reaches  its  maxi- 
iiiiiiii.  iiioie  of  the  l.o<’|’  (  in rent  Mow  lakes  the  southerly 
route,  and  a  COM  ring  is  eventually  pinched  olf. 

The  interesting  point  indicated  by  the  data  pre.sented  here  is 
the  closed  anticvcionic  How  within  the  Loop  Current  immedi- 


Fig  1.5.  Time  history  of  the  rotational  frequency  within  the  Loop 
Current. 


atcly  after  ring  3.378  was  pinehed  o(T.  Also,  Ihe  trajectory  of 
drifter  3354  reveals  that  this  rotational  feature  initially  existed 
o(T  the  coast  of  Cuba.  It  has  been  shown  in  numerical  studies 
that  the  flow  of  the  Loop  Current  can  itself  create  an  anti- 
cyclonic  flow  field  off  the  northwestern  coast  of  Cuba 
IThompsnit,  19861.  The  northward  flow  west  of  the  tip  of 
Cuba  along  with  the  southeasterly  flow  along  the  north  cen¬ 
tral  coast  of  Cuba  will  obviously  produce  negative  vorticity  off 
Cuba's  northwestern  shore.  Such  an  anticyclonic  flow  has 
liceii  documented  by  a  nuntber  of  hydrographic  surveys  [e.g., 
Nouliii  ami  McLellan.  1967;  Cochrane.  1972].  Using  hydro- 
graphic  data.  Hofmann  and  Worley  [1986]  estimated  the  depth 


Fig.  1 7.  'i  ime  hislories  of  Ihe  filtered  observed  speed  components 
for  ring  .1.178  The  light  curve  is  Ihe  easi/west  speed,  and  the  darker 
curve  IS  the  norih/suuth  speed. 
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Till  1  in  nAvT  (  I'lS'l) 

I  ii;  IS  Tiine  hisiiirii-s  of  llie  s«iil  speed  soiiipcitenls  .iboul  the 
seiner  of  rol.ilioii  for  riiii:  TITS  I  he  li^lil  siirse  is  the  east  \scsl  speed, 
and  the  duiker  eiirse  is  the  north  soiilh  speed 

Ilf  this  C'libiin  edd\  In  be  80f)  ni  with  <i  current  ni.ignitude  of 
Id -25  cm  s.  in  good  iigreement  with  the  data  from  this  study 
tl  lutire  12)  Also,  data  slu'w  that  the  temperattirc  and  salinity 
characteristics  of  COM  rings  itnd  the  Cubiin  eddy  are  practi- 
calK  dentical  \ Mrl.Min.  IhfiO.  IMioi.  19821 

I  he  Cuban  eddy  tippcars  to  be  the  How  held  in  which  drif¬ 
ter  ^ '54  initiitlly  hectime  entrained,  fhe  drifter  trajectory  indi- 
s.ites  that  this  Ih'w  held  moved  slightly  northwest  svhile  ring 
"78  was  still  connected  to  the  Loop  Current.  Hut  tt  was  only 
shortiv  after  ring  '.'78  pinched  oil' that  this  anticyclone  pushed 
more  iioithweslerly.  to  iihotil  26  5  N  t-igures  10  and  11  mdi- 
c.iie  III, It  .ilier  this  final  jiiish  (late  August  1985),  the  rotational 
lisdd  iTocame  an  iiilcgial  coin|ionenl  of  the  Loop  Current. 
I  hus  the  data  in  this  case  imply  that  the  kernels  of  Loop 
Cm  lent  rings  may  come  fiom  waters  off  the  northwest  coast 
of  Cuh;i. 


I  le  I')  Time  histories  of  the  east  west  (  Iglit  (.iirsc)  and  north 
soiiih  Id.irkcr  curve)  distances  of  drifter  "7S  (roni  the  center  of  rota¬ 
tion  of  ring  .1.T7X 


Julian  Days  (1985) 

Tig  20  Time  hi.slory  of  the  rotational  frequency  of  ring  3378 

An  example  of  a  modeled  ring  pinch-ofT  in  the  COM  is 
shown  in  I'igures  22  and  23  [from  Walkrali,  1986].  On  day 
63  (Figure  22)  the  model  shows  a  configuration  similar  to  that 
when  drifter  3354  was  seeded.  The  ring  is  well  defined  as  the 
northern  extension  of  the  Loop  Current,  but  it  htis  not  com¬ 
pletely  pinched  oil'  Note  that  there  docs  exist  a  hint  of  anti- 
cyclonic  motion  oir  the  northwestern  coast  of  Cuba.  By  day 
153  (3  months  later)  the  model  shows  the  ring  free  of  the  Loop 
Current.  Moreover,  the  anticyclone  olT  Cuba  is  larger  and 
much  belter  defined.  However,  the  model  shows  the  northern 
edge  of  the  Loop  Current  extending  to  only  26  N.  whereas 
.XBT  and  drifter  data  presented  here  show  an  extension  to 
27  N  in  a  similtir  time  period  (3  months). 

It  is  not  until  after  9  months  before  the  Loop  Current  in 
this  simulation  pushes  northward  again  in  the  process  of 


Julian  Days  (1935) 

l  ig.  21  l  ime  Histories  of  the  eccentricity  (asterisks)  and  major 
,ixis  oricni.iiioii  (dashed  lines)  for  ring  3378  For  the  flow  field  orien- 
i.ilion.  north  is  in  the  positive  e  direction,  and  east  is  in  the  positive 
nine  dircclion. 
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|iii>i.liint:  oil  .iiutllicr  liiip  (rigiire  23),  Uy  (his  lime,  (he  .mli- 
LVelonc  lli.il  h.ul  been  sniin  up  by  (he  Loop  Cm  rent  h.is  lost 
its  ideniily  .I'.iv  '<'>3  fijMire  ’’3)  Thus  (lie  process  ;is  dcler- 
1  '1.1  11  '  •:  "  1  '  oiv  ,il'<;r.l'er  )3S4  i,s  (|i|i<'iciil  from  :his 

,',ii  .cu’ ..  V.  ...il.'l.'o.i  l.i  two  in.ijivi  vv.iys.  I  irst.  the  Loop  Cur- 
•  '  .cd  :  '  I  111!  .'i.iilivvi..i!  to  27  N  in  2  5  months 

I  .if  1(1  rvoiiths  foi  i.fch.iuodcl 

.  .  .  1.  'Cl"'  I'u'.  ..ip'ii  .•\iCiisi''n  was  .lecompaiiied  by 

. ,'iv  ,!  circula  .  voiio  v  orii'iiially  geitcr.ilcd  olf  tlie  coast 

I  V  >ii',.  ili.v. :vci  1.'  olhci  simulations  with  no  deep  in-llovv 
tluouith  the  'itis.ii.in  Straits,  signilicanl  Loop  Current  pene- 
iitiiions  c.iii  ormr  m  40  days  Clearly,  this  is  a  phenomena 
that  rci|uirc.s  further  study 

/.i..'/i  Ctincn!  I'lrsK.v  Rimi  Kiiiciiuitic^ 

Drifters  33.34  and  3378  were  both  in  the  COM  during 
.lulian  davs  IS9  204.  1985.  However,  the  initial  variations  of 
the  kinematic  panimeters  of  ring  3378  (Figures  17-21)  all  indi- 
s.iie  .111  .idiiisimeiil  peiiod  from  days  199  220  during  which 
lime  the  vhiller  (with  its  2()0-m  weighted  line)  found  an  appro- 
pi  itile  orbit  within  the  ring  Thus  we  compare  the  ring  and 
I  0('p  Cm  rent  kineniiilics  from  dtiy  225  to  day  264.  At  the 
beginning  of  Ihi,".  lime  period  the  Loop  Current  was  rotating 
with  I  pernai  of  .iboiil  !0  ^  days,  while  ring  3378  had  a  period 
of  --8  2  vif.y;  '  he  I'rieiitations  of  the  ellipses  of  rotation  (l-ig- 
•ari-  '<1  .Old  21)  were  (|uile  stmil.ir  igenerally  cast  west),  but  the 
L.i’ip  Current  was  sllghtlv  more  elliptical  (e  -  1.5)  than  ring 
377.S  !,■  -  I  2'^' 

I  he  swrl  speeds  tiiu'  (he  radii  of  revolution  for  the  drifters 
III  ri'ili  aiihcvcliii’es  weie  quite  similar  ft'r  vlays  223- 264’  80 
cm  I.’i  ihe  swn!  speed  .  .md  '  91!  km  for  the  i.idii  relative  to 
die  r  .  ’  sei'ter  Inlereitingly,  both  aiitii  VLlones  espeiiented 
'oiieor  I  (.lioCs  of  I'otaiioii  lui  sin.illei  distances  from  the  center 

..I  I  U  i:ii>:i. 

n  "ii'M.  •!,  V  Hchoi’  iiiitl  l/ur  ri  Rim/  Ih-ltichiwni 

We  iif'v  itnn  to  the  kinematics  nf  the  Loop  Cmreiit  at  its 
ail's'  n  utlic'lv  extension  (Julian  davs  230-264)  and  the  kin- 
■l;■.■;ll  of  ring  3378  soon  tiller  it  broke  free  from  the  Loop 
(  III  rent  (Julian  days  225  235)  Kccoing  in  mind  that  these 
d.iM  lepicsent  two  dilVerent  anticyclonic  phenomena,  how 
typical  arc  the  GOM  rings  before  and  after  (hey  break  olF 
from  the  Loop  Current? 

The  similarities  between  the  two  sets  of  kinematics  arc  con¬ 
siderable.  First,  the  magnitudes  of  the  swirl  speeds  and  of  the 
radii  arc  practiciilly  identical.  Moreover,  their  elliptical  orien¬ 
tations  arc  both  northvvcst/southcasl.  The  Loop  Current  was 
shgh'.lv  more  elliptical  than  ring  3378.  but  one  might  expect 
-iieli  .1  dillercncc  in  ellipticiiy  seeing  that  an  anticyclone  still 
.'ii.ivhed  to  the  Loop  Current  is  a  "captmed"  phenomena, 
while  .1  detached  iing  is  an  isolated  vortex.  The  similarities  of 
the  .iniiey clones  plus  their  close  periods  of  rotation  imply  that 
tl’.c  basic  kinematic  characteristics  of  GOM  rings  can  be  es¬ 
tablished  as  the  Loop  Current  pushes  northwestward  olF  the 
slu'rc  of  Cuba. 

Siininuirv 

rile  scenario  implied  here  is  one  in  which  the  essence  of  the 
ling  IS  esttiblishcd  oil' the  northwest  cotist  of  Cuba,  develops  in 
deeper  waters  oil'  ihe  coastline,  and  then  takes  on  its  final 
characteristics  shortly  after  the  previous  ring  coiiipletely  de¬ 
taches  from  the  I  oop  Current. 


The  implication  of  such  a  mechanism  is  that  of  time  scales 
One  is  not  required  to  wail  until  instabilities  in  the  (low  field 
grow  larcc  enough  to  produce  a  closed  rotational  feature.  The 
data  heie  show  that  a  iica  Loop.  Cnirmt  rotational  leattire 
can  be  well  established  even  be.’'oic  the  previous  ring  has  to¬ 
tally  pinched  oil'.  It  was  cstiiblidicd  that  line  3378  was  de¬ 
tached  from  the  Loop  Curici'.:  by  ii:id-Jii!y  l';85.  Vet  It  was 
seen  that  the  ne.v  rotatioiml  stiucture  had  moved  to  latitude 
26.5'N  by  the  beginning  of  September  1985  Thii.s  the  Loop 
Current  was  rcconligured  for  another  ring  separation  only  1.5 
months  after  the  previous  ring  sc|niratioii.  Elliot  [1982]  docu- 
mefitcd  three  ring  separations  in  a  12-nionth  period  using 
hydrograpliic  data.  The  presence  of  Ghost  Eddy  along  with 
ring  3378  implies  a  similar  separation  rate.  Model  studies 
[c.g.,  Hiirlburl  ami  Thompson,  1980;  Thompson,  1986;  Wall- 
craft,  1986]  report  separation  rates  ranging  from  4  to  18 
months,  depending  upon  the  type  of  instability  involved. 
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Observed  and  Simulated  Kinematic  Properties  of  Loop  Current  Rings 

A.  D.  Kirwan.  J.  K.  Lewis, ^  A.  W.  Indest,'*  P.  Reinersman,‘  and  I.  Quintero* 

Two  rings,  shed  by  the  Loop  Current  in  1980  and  1982.  were  observed  for  several  months  by 
satellite-tracked  drifters  to  migrate  across  the  Gulf  of  Mexico.  The  drifter  path  data  have  been  inverted 
to  obtain  estimates  of  the  paths  of  the  centers  of  the  two  nngs.  nng  shape,  and  the  swirl  velocities.  Three 
drifters  were  deployed  in  the  1980  nng.  and  the  analysis  of  that  data  set  establishes  the  vanability  of  the 
above  kinematic  estimates  for  one  nng.  A  companson  of  the  analysis  of  data  from  both  rings  provides 
some  idea  on  mter-ring  vanability.  Both  nngs  impacted  the  Mexican  continental  slope  at  about  22.8°N. 

95  5'W  After  a  brief  adjustment  penod.  both  rings  reestablished  and  maintained  a  vortex  character  for 
several  months  in  the  slope  region  while  migrating  slowly  to  the  north.  The  paths  of  the  centers  of  the 
two  rings  along  the  slope  are  virtually  identical.  The  same  analysis  routine  was  applied  to  some  simulat¬ 
ed  dnfter  data  obtained  from  the  Hurlburt  and  Thompson  (1980)  Gulf  of  Mexico  primitive  equation 
model.  In  the  midgulf,  the  agreement  between  the  observed  nngs  and  the  simulated  ring  is  good, 
although  the  former  showed  stronger  interaction  with  the  continental  slope  topography  and/or  circu¬ 


lation  than  was  seen  in  the  latter.  Along  the  slope, 
extraordinary  agreement  with  the  observations. 


1.  Introduction 

The  circulation  of  the  Gulf  of  Mexico  has  been  the  subject 
of  a  surprising  amount  of  theoretical  and  observational  stud¬ 
ies.  From  the  turn  of  the  century  [^Sweitzer,  1898]  until  1973 
[,4ustin.  1955;  Nowlin  and  McClellan,  1967;  Nowlin,  1972],  the 
emphasis  was  on  analyses  of  hydrographic  data.  These  studies 
established  the  presence  of  poorly  defined,  but  quasi¬ 
permanent,  anticyclonic  signatures  in  the  central  part  of  the 
gulf. 

Since  1973,  the  emphasis  has  shifted  from  qualitative  de¬ 
scriptions  of  the  hydrography  to  attempts  at  quantifying  dy¬ 
namical  mechanisms  associated  with  the  anticyclonic  struc¬ 
tures.  Siurges  and  Blaha  [1976]  and  Blaha  and  Stnrges  [1978] 
proposed  that  the  western  gulf  circulation  was  the  result  of  the 
balance  of  wind  stress  curl  and  planetary  vorticity,  i.e.,  a  mini¬ 
western  boundary  current  effect  in  the  Gulf  of  Mexico.  Elliot 
[1979,  1982]  disputed  this,  noting  that  this  balance  was  not 
generally  achieved  if  synoptic  wind  data  on  a  fine  scale  were 
used  His  work  also  gave  more  credence  to  the  earlier  idea  of 
Ichiye  [1962]  that  the  anticyclonic  features  were  in  fact  rings 
that  had  been  shed  by  the  Loop  Current  in  the  eastern  gulf. 
Kirwan  et  al.  [1984a]  substantiated  this  hypothesis  by  track¬ 
ing  the  migration  of  a  ring  shed  by  the  Loop  Current  across 
the  gulf  to  the  continental  shelf  off  of  Mexico. 

In  a  companion  paper,  Kirwan  et  al.  [19846]  assessed  the 
translation  velocities,  local  vorticity,  deformation  rates,  and 
shape  of  the  ring  as  it  propagated  across  the  gulf.  The  analyses 
of  the  three  drifters  that  were  in  the  ring  simultaneously 
showed  generally  good  agreement  of  these  properties.  The  re¬ 
sults  also  are  in  general  qualitative  agreement  with  the  nu- 
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the  model  ring  kinematic  charactenstics  were  m 


merical  simulations  from  the  general  circulation  model  (GCM) 
of  Hurlburt  and  Thompson  [1980],  hereinafter  referred  to  as 
HT. 

Is  the  qualitative  agreement  between  GCM  calculations  and 
one  particular  ring  fortuitous  or  the  result  of  deterministic 
dynamics?  We  examine  this  question  with  observations  from 
two  different  rings  along  with  simulated  data  from  the  HT 
model.  Specifically,  we  compare  in  detail  kinematic  properties 
such  as  the  movement  of  ring  centers,  ring  translation  and 
swirl  velocities,  and  ring  shapes  from  1980-1981,  1982-1983, 
and  simulations  from  the  HT  GCM. 

The  analysis  routine  used  to  determine  these  kinematic 
properties  is  an  improved  version  of  that  reported  by  Kirwan 
et  al.  [19846].  The  routine  inverts  Lagrangian  path  data  to 
obtain  the  desired  kinematic  properties.  Unlike  the  earlier 
study,  the  present  routine  emphasizes  ring  shape  and  trans¬ 
lation.  Also,  considerably  better  time  resolution  on  parameter 
estimates  can  be  made  with  the  new  routine.  Some  details  are 
given  below  and  in  the  appendix.  The  internal  consistency  of 
the  analysis  routine  is  being  assessed  separately  with  simulat¬ 
ed  data  from  the  HT  GCM. 

The  following  observed  data  were  utilized  in  this  study. 
These  come  from  two  rings  which  occurred  in  1980-1981  and 
1982-1983.  We  have  reported  previously  on  the  former  [see 
Kirwan  et  al.,  1984a,  6]  using  the  more  primitive  analysis 
routine.  That  nng  had  three  drifters  in  it  simultaneously,  and 
so  our  re-analysis  provides  some  indication  of  intraring  varia¬ 
bility  along  with  some  indication  of  the  consistency  of  the 
analysis  routine.  Only  a  descriptive  analysis  of  the  latter  ring, 
which  had  only  one  drifter,  has  been  reported  before  [Lewis 
and  Kirwan,  1985].  Comparison  of  these  two  rings  provides 
some  measure  of  inter-ring  variability. 

These  observed  data  are  compared  with  simulated  La¬ 
grangian  data  from  the  HT  GCM.  This  model  was  a  two-layer 
(lower  layer  active)  nonlinear  primitive  equation  model  of  the 
Gulf  of  Mexico  with  realistic  bottom  topography.  The  grid 
spacing  for  the  calculations  used  here  was  10  km.  The  simulat¬ 
ed  path  was  obtained  by  advecting  a  parcel  at  each  model 
time  step  (4  hours)  to  a  new  position  given  by  the  product  of 
the  velocity  vector  at  the  current  position  of  the  parcel  and  the 
time  step.  Since  the  position  of  the  parcel  was  rarely  at  a  grid 
point,  it  was  usually  necessary  to  interpolate  t!ie  velocity 
vector  from  surrounding  grid  points.  These  calculations  were 
kindly  performed  by  A.  Wallcraft. 
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Both  the  observations  and  the  simulations  track  appropri¬ 
ate  rings  for  many  months.  During  this  time,  the  real  and 
simulated  rings  move  from  the  midgulf  region  with  fairly  flat 
topography  in  excess  of  3000  m  to  the  continental  slope  off  of 
Mexico  This  area  of  interest  is  depicted  in  Figure  1.  In  the 
midgulf  region,  ring  translations  are  governed  essentially  by 
Rossby  wave  dynamics,  and  so  one  expects  them  to  move 
west-southwest  But  along  the  slope,  interactions  of  the  ring 
with  bottom  topography  become  important.  Since  the  dynam¬ 
ics  in  these  two  geographic  regions  are  so  different,  we  have 
performed  separate  analyses  for  the  rings  in  the  two  regions. 


Fig  I  Depiction  of  the  area  in  which  the  observed  and  simulated 
Loop  Current  rings  are  studied 


2.  Analysis  Procedure 

The  observations  used  in  this  study  are  the  positions  of  the 
Polar  Research  Laboratory  drifters  with  Argos  identifications 
1598,  1599,  1600,  and  3374.  Kirwan  et  al.  [1984a,  b]  reported 
on  the  analysis  of  1598,  1599,  and  1600,  which  were  deployed 
simultaneously  in  November  1980.  Position  data  from  these 
three  drifters  are  analyzed  to  establish  the  variability  within  a 
ring.  Intraring  variability  is  assessed  by  comparing  the  analy¬ 
sis  performed  on  the  three  drifters  in  the  1980-1981  ring  with 
that  of  drifter  3374,  which  was  deployed  in  the  1982-1983 
ring.  Each  drifter  was  the  standard  Polar  Research  Labora¬ 
tory  buoy  drogued  by  a  200-m  thermistor  line.  No  thermistor 
data  are  available  from  any  of  these  units. 

The  analysis  procedure  utilizes  a  parametric  kinematic 
model  discussed  by  Kirwan  et  al.  [1984b]  for  inferring  the 
translation  and  swirl  velocities  as  well  as  the  vorticity,  hori¬ 
zontal  divergence,  and  deformation  rates.  This  model  was  first 
proposed  by  Okubo  [1970]  to  describe  small-scale  motion 
about  flow  singularities.  Kirwan  [1984]  showed  that  same 
model  can  be  applied  to  larger  scales  if  the  model  parameters 
are  constant  along  a  parcel  path. 

The  model  assumes  that  the  swirl  velocity  is  given  as  a 
linear  function  of  the  distance  from  a  local  flow  center.  Here 
the  local  parcel  flow  center  will  be  referred  to  as  the  ‘‘ring 
center”  or  some  times  "orbit-determined  center”.  This  termin¬ 
ology  is  dictated  only  by  the  nature  of  the  available  data. 


Flierl  [1981]  has  pointed  out  that  this  center  may  not  coincide 
with  an  "Eulerian”  center.  (We  interpret  this  as  a  center  of 
mass.)  Smith  and  Reid  [1982]  have  developed  a  general  for¬ 
malism  for  studying  the  centers  of  the  distribution  on  any 
property,  for  example,  mass,  potential  energy,  kinetic  energy 
or  enstrophy.  In  general,  these  centers  do  not  coincide,  nor  do 
their  propagation  velocities.  Establishing  the  relations  be¬ 
tween  these  various  centers  for  Gulf  of  Mexico  rings  requires 
much  more  daf*  than  is  now  available. 

For  the  present  model  the  total  velocity  (u,  u)  of  a  drifter  is 
composed  of  the  translation  velocity  [Uf,  of  the  ring 
center  and  the  swirl  velocity  (Ms,  Uj); 

M  =  t/r  +  uj 

(1) 

u  = 

Us  =  {d  +  a)x/2  +  (b  -  c)y/2 

(2) 

Vs  =  (b  +  c)x/2  +  {d-  a)y/2 

Here  x  and  y  are  the  instantaneous  coordinates,  relative  to  the 
ring  center  of  a  particular  parcel  (drifter).  If  (2)  is  regarded  as  a 
Taylor  expansion  for  the  velocity  field,  then  the  parameters  a 
and  b  can  be  interpreted  as  deformation  rates  and  c  and  d  as 
the  vertical  vorticity  and  horizontal  divergence.  It  is  stressed 
that  these  parameters  are  evaluated  along  individual  paths. 
Eulerian  estimates  of  deformation  rates,  vorticity,  and  diver¬ 
gence  likely  would  differ  from  these  estimates.  Consequently,  it 
is  better  to  interpret  a,  b,  and  c  as  shape  functions  which 
determine  the  orientation  and  ellipticity  of  a  particular  orbit 
in  a  ring.  See  Kirwan  et  al.  [1984b]  for  details. 

The  general  solution  to  (1)  and  (2)  is  easily  obtained 
lOkubo,  1970].  The  procedure  here  is  to  apply  this  solution  to 
every  time  interval  between  fixes  (in  practice,  interpolated 
positions).  Thus  for  the  interval  t^  ^  t  ^  t^^.^,  one  obtains 

+  (dt  +  o»)V2  +  (h»  -  Ct)y»/2  (3) 

I'l  =  ^Tk  +  +  Ck)xJ2  +  (d,  -  Ot)yj/2  (4) 

where 

x»  =  {exp  [ru(t-l»)][2rt(yV  +  V-^+flE)  li(^-Ct)] 
-exp  [ra(f  -  f»)]  [  Aftfaj  -  v/o»  ^  +  b*  *  -  c*  ^)  -F  Ttfbj-Ci)]} 

-  27a7+V^^  (5) 

y*  =  {exp[ru(i-g][2f»(b*  +  ct)-l-  Waj^ +  bt^ -Ct^-a*)] 

+exp  ['■2*(f  -  ft)]  [  -  Ok^  +  bk^-Ck^+ ''»)] } 

-  2^a^^  +  bj^  -  (6) 

In  (5)  and  (6),  (2fj,  Y^)  are  the  coordinates  of  the  parcel  in 
question  at  time  t  =  tj  relative  to  the  ring  center.  Also, 

r.*  =  (d*-F^47Tb7^)/2  (7) 

rik  =  (<i»  -  s/ok^  +  bk^  -Ck^)/2  (8) 

are  the  eigenvalues  of  the  matrix 

(dt  +  Oj)  (bt-Cj)"! 

ibk  +  Ct)  (<it-ajJ 

Note  that  these  eigenvalues  are  complex  conjugates  when 
Cj^  >  +  bj*.  This,  naturally,  produces  real  periodic  solu- 
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lions  which  are  a  characteristic  of  swirl  velocities  of  drifters 
trapped  in  a  ring. 

To  apply  this  model  to  drifter  data,  diflferentiate  (3)  and  (4) 
three  times  with  respect  to  time.  One  side  of  these  equations 
contains  the  unknown  variable  set  Zj  =  {a^,  b^,  d^,  X,,,  li). 

The  other  side  of  these  equations  can  be  estimated  from  path 
data.  These  six  equations  do  not  contain  the  translation  veloc¬ 
ity.  since  it  is  assumed  a  priori  to  be  constant  for  a  time 
interval.  (But,  of  course,  its  value  may  change  from  one  time 
interval  to  the  next.)  These  six  equations  are  inverted  analyti¬ 
cally  for  the  six  elements  of  Z,,.  The  translation  velocities  can 
be  recovered  from  (3)  and  (4),  smce  all  other  quantities  are 
now  known.  Alternatively,  the  translation  velocities  can  be 
obtained  by  differentiating  the  path  of  the  ring  center.  Both 
procedures  provide  comparable  estimates.  Here  we  use  the 
latter  procedure.  With  all  elements  of  now  determined,  one 
then  moves  to  the  next  time  interval  and  repeats  the  calcula¬ 
tions,  thus  obtaining  a  time  senes  of  Z^.  Details  are  given  in 
the  appendix. 

We  have  expended  a  great  deal  of  effort  in  learning  how  to 
make  reliable  time  derivative  estimates  of  the  path  data.  For 
the  data  analyzed  here,  the  raw  position  files  have  been  edited, 
interpolated  to  equally  spaced  time  intervals,  and  low-pass 
filtered  (half  power  point  of  1(X)  hours).  All  derivative  esti¬ 
mates  were  obtained  by  fourth-order  accurate,  centered  differ¬ 
ence  schemes.  Each  derivative  file  was  edited  and  low-pass 
filtered.  The  calculations,  as  described  in  the  appendix  for  Z^, 
were  performed  with  the  latter  derivative  files.  The  calculated 
Z»,  in  turn,  were  edited  and  low-pass  filtered. 

The  analysis  procedure  just  outlined  is  an  improvement 
over  that  used  by  Kirwan  et  al.  [19846].  Their  method  re¬ 
quired  an  a  priori  assignment  of  a  solution  form  to  (4)  and  (S); 
i.e..  (7)  and  (8)  were  assumed  to  be  complex  conjugates,  and 
values  of  Z^  were  assumed  to  be  constant  over  one  ring  revo¬ 
lution  (approximately  15  days).  The  procedure  used  here  does 
not  require  an  a  priori  assignment  of  the  solution  form;  it  lets 
the  data  determine  this.  Thus  it  should  have  wider  utility.  In 
addition,  the  Z^  are  now  required  to  be  constant  over  seven 
time  intervals  (42  hours  in  the  present  case)  rather  than  2 
weeks. 

As  was  indicated  above,  a,  b,  and  c  are  not  especially  useful 
for  comparison  with  Eulerian  data.  For  rings  a  more  appro¬ 
priate  description  would  be  the  elliptical  structure  of  the  orbit 
traversed  by  a  particular  drifter.  As  was  shown  by  Kirwan  et 
at.  [19846],  the  characteristic  equation  is 

(Ct  +  6*)xt^  -t-  (Cj  -  h^y^^  -  2a»Xiyj 

=  Ljexp[-d*(f-fi)]  (9) 

Here  is  the  angular  momentum  per  unit  mass  relative  to 
the  ring  center.  This  can  b:  calculated  directly  (see  equation 
(A  18)).  In  (9)  the  argument  in  the  exponential  tenn  can  be 
made  zero  by  making  the  evzJuation  at  the  beginning  of  the 
time  interval  t  =  t^.  Note  that  a  requirement  for  (9)  to  describe 
an  ellipse  is  that  -I-  6*^. 

From  analytic  geometry  it  is  known  that  for  anticyclonic 
motion,  c^  <  0,  the  major  axis  of  the  ellipse  makes  an  angle 
with  east  of 

at*|tan‘‘ (-aj/6j)  (10) 

In  evaluating  (10),  care  must  be  exercised  in  the  quadrant 
assignment.  Also,  the  semilength  of  the  major  axis  is 


—  W  ^Jlfik  ~  ( 1 U 

where 

hk^=a,^+b,^ 

The  paradigm  just  outlined  for  analysis  of  Lagrangian  data 
puts  critical  emphasis  on  time  derivatives  of  path  data.  Incon¬ 
sequential  errors  in  the  path  data  can  become  consequential  in 
the  derivative  estimates.  Moreover,  the  nonlinear  algebraic 
combinations  of  these  derivatives  in  the  paradigm  may  mag¬ 
nify  further  the  impact  of  these  errors  on  the  inversion.  The 
concern  is  that  two  paths  which  contain  the  same  kinematic 
information  but  differ  by  small  random  and/or  round  off 
errors,  will  produce  vastly  different  kinematic  estimates  when 
run  through  the  paradigm.  The  internal  consistency  of  esti¬ 
mates  of  ring  kinematics  is  presently  being  examined  with 
simulated  data  from  the  HT  GCM. 

3.  Midgulf 

3.1.  Observations 

Here  the  results  of  the  intraring  and  inter-ring  comparisons 
are  presented  for  the  midgulf  region.  This  includes  ring  kin¬ 
ematics  as  inferred  from  the  three  drifters  in  the  1980-1981 
ring  and  the  single  drifter  in  the  1982-1983  ring.  While  in  the 
midgulf  region,  both  rings  began  to  interact  with  the  conti¬ 
nental  slope  topography  and/or  circulation.  This  period  is  iso¬ 
lated  in  the  analysis.  The  same  analysis  routine  is  then  applied 
to  simulated  midgulf  ring  data  from  the  HT  GCM. 

The  1980-1981  ring.  The  data  for  this  portion  of  the  study 
come  from  three  drifters  (1598,  1599,  and  1600)  which  were 
seeded  in  the  198(V-1981  ring.  Kirwan  et  al.  [1984o,  6]  have 
already  reported  on  this.  The  data  have  been  reanalyzed  using 
the  general  algorithm  discussed  in  the  preceding  section  and 
in  the  appendix. 

Figure  2  shows  the  paths  of  these  three  drifters  in  the 
midgulf  region  along  with  the  positions  of  the  centers  and  the 
orbit  ellipses.  Arrows  on  the  drifter  and  center  paths  mark 
10-day  intervals,  while  the  ellipses  are  presented  at  15-day 
intervals.  The  paths  for  each  drifter  are  for  a  common  time 
interval  (day  340,  1980,  to  day  118,  1981).  This  figure  es¬ 
tablishes  three  points.  First,  all  three  make  the  same  swath 
through  this  portion  of  the  gulf,  suggesting  that  they  generally 
followed  the  ring.  The  overall  path  characteristics  of  drifters 
1598  and  1600  (Figures  2a  and  2c)  are  remarkably  similar, 
suggesting  that  they  were  in  very  similar  orbits.  Second,  the 
paths  for  the  center  of  the  ring,  as  determined  independently 
from  the  three  drifters,  is  in  good  agreement  as  far  as  about 
94°W  (the  first  2  months  of  the  record).  Discounting  the  very 
beginning  where  there  are  numerical  problems,  the  typical  dif¬ 
ference  in  the  location  of  the  center  at  any  one  time  is  30  km, 
about  20%  of  the  diameter  of  the  ring  as  determined  from 
satellite  imagery  [Kirwon  et  al.,  19846].  This  is  about  10  km 
larger  than  the  rms  di.sp' icement  between  the  maximun  sur¬ 
face  and  interface  displacement  anomalies  of  the  HT  GCM. 
The  30-km  variability  is  considerably  larger  than  that  report¬ 
ed  by  Hooker  and  Olson  [1984]  using  a  variant  of  the  tech¬ 
nique  employed  here  under  ideal  conditions.  In  view  of  our 
results  along  the  continental  slope,  described  below,  it  is  not 
clear  whether  the  large  variability  is  due  to  data,  technique,  or 
rapid  evolution  of  the  ring. 

Finally,  the  ellipses  also  show  generally  good  agreement 
east  of  94°W.  Note  the  tendency  to  develop  a  northwest  elon- 
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gation  as  the  ring  approaches  94^W  West  of  94°W,  all  three 
orbits  suggest  that  the  ring  begins  to  interact  with  the  conti¬ 
nental  slope  circulation  and/or  topography.  This  is  discussed 
later. 


98»  97“  96*  95°  94°  93'  92»  gi"  90" 


LONGITUDE 
Fig.  2a 


Figure  3  shows  the  time  series  for  the  swirl  and  translation 
velocities.  All  three  records  indicate  swirl  velocities  in  excess  of 
50  cm  s  throughout  most  of  the  record.  The  dominant  peric  d 
IS  about  14  days  with  a  variation  of  about  2  days  between  the 
records.  This  is  recognized  as  the  ring’s  rotation  period.  How¬ 
ever.  from  about  day  30  to  60.  1981.  all  three  drifters  show  a 
smaller-period  component  (about  5  days)  in  the  time  series 
and  a  decrease  in  swirl  amplitude.  As  will  be  discussed  below, 
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Fig.  2.  Drifter  paths,  inicrred  center  path  (red)  and  orbit  ellipse 
(blue)  for  the  three  dnfters  in  the  1980-1981  ring;  (a)  1598,  (6)  1599, 
and  (c)  1600.  Arrows  on  the  paths  are  at  10-day  intervals,  while  the 
ellipses  are  shown  at  1 5-day  intervals.  Here  S  and  F  refer  to  start  and 
finish. 
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Fig.  3.  Swirl  (left  scale)  and  translation  velocities  (nght  scale)  in¬ 
ferred  from  the  three  dnfters  in  the  1980-1981  nng.  (a)  1598,  (6)  1599, 
and  (c)  16<X).  The  east  and  north  components  of  the  translation  veloc¬ 
ity  are  depicted  in  red  and  blue,  respectively. 
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It  IS  speculated  that  during  this  period  the  drifters  were  re¬ 
sponding  to  effects  produced  by  the  continental  slope  and,  or 
circulation  further  to  the  west.  All  the  translation  velocities 
show  a  consistent  westward  component  of  about  4  cm,s  up 
until  about  day  30.  1981.  For  the  reason  just  given,  it  is  not 
clear  that  the  analysis  applies  to  this  ring  for  the  period  from 
day  30.  1981.  to  day  60,  1981  After  day  60,  1981,  the  ring 
centers  appear  to  be  stationary. 

The  1982  1983  nnq  and  comparison  with  the  1980-1981 
ring.  The  same  characteristics  are  seen  in  the  analysis  for  the 
1982-1983  ring.  Figure  4  shows  the  drifter  path,  inferred 
center  path,  and  the  orbit  ellipses  for  drifter  3374.  As  in  Figure 
2.  there  is  significant  distortion  of  the  orbits  west  of  94‘’W. 
Unlike  the  earlier  ring,  no  significant  northwest  ellipticity  is 
developed  as  the  ring  approaches  94'"W.  Comp.irison  with 
Figure  2  indicates  that  the  1982-1983  ring  center  followed 
nearly  the  same  path  as  that  of  the  1980-1981  ring. 


Fig.  4.  Drifter  path,  inferred  center  path  (red),  and  orbit  ellipse 
(blue)  for  the  drifter  in  the  1982-1983  ring.  Arrows  on  the  paths  are  at 
10-day  intervals,  while  the  ellipses  are  shown  at  IS-day  intervals. 
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Fig.  5.  Swirl  (left  scale)  ana  translation  velocities  (right  scale)  in¬ 
ferred  from  the  three  drifters  in  the  1982-1983  nng.  The  east  and 
north  components  of  the  translation  veloaty  are  depicted  in  red  and 
blue,  res.,jctively. 

Figure  5  shows  the  translation  and  swirl  velocities  for  the 
1982-1983  ring.  The  swirls  are  slightly  less  than  50  cm/s, 
which  could  reflect  a  smaller  orbit  of  the  drifter.  The  dominant 
period  IS  still  about  14  days.  Note  that  around  day  1,  1983  the 


swirl  velocities  decrease  significantly  and  show  some  evidence 
of  higher-frequency  components.  This  is  the  penod  when  the 
drifter  moves  west  of  94°W  and  executes  the  deformed, 
northwest-oriented  loop.  As  with  the  earlier  ring,  this  is  inter- 
f  reted  as  evidence  for  topographic  and/or  slope  circulation 
interactions. 

Figure  2  indicates  that  west  of  94°W,  the  three  drifters  in 
the  1980-1981  ring  show  significant  divergence  in  the  calcula¬ 
tion  of  the  inferred  center  Furthermore,  drifters  from  both 
rings  showed  anomalous  path  characteristics  in  this  region.  In 
particular,  1598  and  16(X)  appeared  to  become  entrained  in 
smaller,  anticyclonic  eddies,  while  1599  and  3374  exhibited 
warped  and  elongated  ellipses.  This  has  been  attributed  to 
interaction  between  the  ring  and  the  continental  slope  and/or 
slope  circulation  further  to  the  west. 

Figure  6  focuses  on  this  time  interval.  Figure  6a  shows  the 
paths  of  all  four  drifters  during  the  period  in  question.  It  is 
clear  from  this  that  all  four  drifters  were  under  the  influence  of 
different  dynamical  processes  here  than  those  they  experienced 
to  the  east. 


Fig.  6<i 


Fig.  6.  Paths  of  drifters  1598,  1599  and  1600  in  the  1980-1981 
nng  and  3374  in  the  1982-1983  nng  for  the  penod  when  they  were 
interacting  with  the  slope  topography  and/or  circulation.  The  arrows 
are  at  10-day  intervals,  (a)  Drifter  paths,  (b)  Inferred  centen  of  flow. 

Figure  6b  compares  the  movement  of  the  inferred  cen.ers. 
In  interpreting  this  figure,  it  should  be  emphasized  that  the 
center  calculation  is  greatly  complicated  by  the  changes  in 
curvature  of  the  paths.  The  paradigm  interprets  this  as  a 
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change  from  an  anticyclonic  to  cyclonic  structure  with  a 
consequent  change  in  the  location  of  the  center.  No  doubt  the 
real  dynamics  during  this  period  are  much  more  complicated 
than  the  rather  simple  scenario  available  from  the  kinematic 
model.  The  geographic  variation  in  the  paths  clearly  is  much 
larger  than  that  seen  east  of  94°W  or.  as  will  be  seen  shortly, 
than  that  found  along  the  slope.  This  suggests  that  each  of 
these  drifters  was  either  briefly  pulled  out  of  the  parent  ring  or 
its  orbits  deformed  through  interaction  with  the  continental 
slope  topography  or  circulation.  This  hypothesis  accounts  for 
the  deformed  orbits  and  the  brief  occurrence  of  high  fre¬ 
quencies  in  the  swirl  velocities. 

3.2.  Simulations 

Attention  is  now  turned  to  the  simulated  data.  There  are 
two  issues  involved  in  the  utilization  and  interpretation  of  this 
data.  First,  how  consistent  is  the  paradigm  in  recovering  ring 
kinematics  which  should  be  independent  of  the  orbit?  This 
issue  is  presently  being  addressed.  Addressed  here  is  the  issue 
of  establishing  how  well  the  GCM  agrees  with  the  observa¬ 
tions.  Speciflcally,  the  center  path,  orbit  ellipses,  and  trans¬ 
lation  velocity  as  inferred  from  the  simulated  data  are  com¬ 
pared  with  the  same  properties  as  determined  from  the  obser¬ 
vations. 
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Figure  7  shows  the  path  of  one  of  several  simulated  drifters 
from  the  HT  GCM  along  with  the  orbit  ellipses.  Figure  7a 
depicts  the  results  using  the  ring  center  as  inferred  from  the 
drifter  orbit,  while  Figure  7b  repeats  the  ellipse  calculation 
using  the  position  of  the  maximum  displacement  anomaly  of 
the  interface.  The  choice  of  the  interface  displacement  as  rep¬ 
resentative  of  the  ring  center  is  arbitrary.  Presently,  we  are 
comparing  the  maximum  interface  displacement,  the  maxi¬ 
mum  surface  elevation  displacement  and  orbit-determined 
centers.  In  general,  all  are  within  30  km  of  each  ether. 

The  most  significant  difference  between  the  orbit- 
determined  centers  and  interface  displacement  anomaly  occurs 
in  the  interval  between  91°W  and  93°W  where  the  orbit- 
inferred  ring  center  shows  a  significant  southern  loop.  This  is 
not  seen  in  the  path  of  the  maximum  interface  displacement 
anomaly.  At  93'“W  the  orbit-inferred  center  executes  an  anti¬ 
cyclonic  loop  that,  again,  is  not  seen  in  the  displacement 
anomaly  path.  Elsewhere,  the  orbit-determined  path  is  consis¬ 
tently  5  to  25  km  south  of  the  displacement  anomaly  path. 

A  comparison  of  Figures  2,  4,  and  7  shows  that  the  simulat¬ 
ed  and  observed  centers  follow  the  same  general  path  through 
this  part  of  the  midgulf.  In  general,  the  simulated  center  path 
is  10  to  35  km  south  of  the  observed  path,  depending  upon 
which  of  the  two  simulated  center  paths  is  used.  There  is 
virtually  no  evidence  of  the  interaction  with  the  continental 
slope  topography  and  circulation  in  the  model  tesults  as  was 
seen  in  the  observations.  None  of  the  simulated  orbits  were 
deformed  or  appeared  to  be  pulled  out  temporarily  from  the 
ring.  Somewhat  surprisingly,  we  have  found  less  consistency  in 
the  ellipse  calculations  for  the  simulations  than  in  the  observa¬ 
tions.  A  number  of  the  ellipses  from  the  simulated  data  are 
quite  deformed.  They  also  show  frequent  reversals  of  the 
major  and  minor  axes.  Overall,  however,  as  the  ring  nears  the 
slope  region,  a  northeast  orientation  appears  to  develop  in  the 
orbit  ellipses. 

Figure  8  shows  the  translation  and  swirl  velocities  for  the 
simulation.  Figure  8u  gives  the  results  using  the  orbit  center, 
and  Figure  8b  uses  the  maximum  displacement  anomaly  as 
the  ring  center.  The  dominant  period  in  the  swirl  velocities  is 
about  24  days,  or  10  days  more  than  was  found  in  the  obser¬ 
vations.  Note  that  the  swirl  velocities,  as  calculated  from 
cither  center,  are  of  comparable  magnitude,  about  30  cm/s. 
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Fig.  7b 

Fig.  7.  Path  and  orbit  ellipses  (blue)  of  a  simulated  drifter  from 
the  HT  GCM.  Arrows  on  the  path  are  at  10-day  intervals,  while 
ellipses  are  shown  at  1 5-day  intervals.  Figure  7a  shows  the  ellipses 
using  the  orbit-inferred  center  as  the  flow  center,  and  Figure  7b  uses 
the  maximum  interface  displacement  anomaly.  Both  centers  are 
shown  in  red. 


Fig.  8fl 

Fig.  8.  Translation  and  swirl  velocities  as  inferred  from  the  HT 
GCM.  The  east  and  north  components  of  the  translation  velocity  are 
depicted  in  red  and  blue,  respectively.  The  scale  for  the  translation 
velocity  is  on  the  nght.  Figure  8a  shows  results  using  the  orbit- 
inferred  center,  while  Figure  8b  {next  page)  shows  results  using  the 
maximum  displacement  anomaly  as  the  flow  center. 
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Fig.  8b 


This  is  somewhat  less  than  the  observed  swirls,  but  this  could 
be  rectified  by  picking  a  simulated  drifter  in  a  greater  orbit. 
There  seems  to  be  more  difference  in  the  translation  velocities. 
Because  its  path  is  straighter,  the  displacement  anomaly  path 
has  less  variance.  The  mean  westward  component  in  both 
cases  of  about  5  cm/s  is  consistent  with  the  observations. 
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Fig.  9.  The  centers  for  the  1980-1981  ring,  the  1982-1983  ring, 
and  the  HT  GCM  simulation.  Orbit  I  (red)  refers  to  the  orbit- 
determined  center,  and  pressure  (blue)  refers  to  the  maximum  dis¬ 
placement  anomaly  path.  Arrows  depict  10-day  intervals. 

Figure  9  compares  the  inferred  center  movement  for  all  four 
drifter.'  along  with  the  inferred  center  of  the  simulated  drifter 
and  the  displacement  anomaly  path  east  of  94°W.  This  is  the 
region  where  there  is  no  apparent  effect  of  slope  topography 
or  circulation  in  the  observations.  Of  course,  details  in  the 
paths  vary,  but  overall  the  agreement  is  quite  good.  The  maxi¬ 
mum  deviation  at  any  one  time  between  any  of  the  paths  is 


100  km.  which  is  only  two-thirds  the  diameter  of  a  typical 
ring.  For  most  of  the  paths,  the  orbit-determined  centers  are 
somewhat  south  of  the  displacement  anomaly  path.  However, 
as  the  paths  near  94°W,  the  observed  centers  move  north  of 
the  displacement  anomaly.  All  the  orbit-determined  center 
paths  show  considerably  more  structure  than  does  that  of  the 
displacement  anomaly. 

4.  Continental  Slope 

It  is  remarkable  that  both  observed  rings  and  the  simulated 
ring  impacted  the  continental  slope  in  the  same  area.  This  is 
seen  in  Figure  10,  which  shows  the  drifter  paths,  inferred 
center  paths,  and  orbit  ellipses  for  all  three  cases.  In  the  case 
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Fig.  10.  Drifter  paths,  inferred  centers  and  orbit  ellipses  for  (a)  the 
1981-1982  ring,  (b)  the  1982-1983  ring  (above);  (c)  the  simulated  ring, 
and  (d)  center  paths  along  the  continental  slope  (next  page).  Arrows 
on  the  paths  depict  10-day  intervals,  while  ihe  ellipses  are  at  15-day 
intervals. 

of  the  simulated  data  (Figure  10c),  the  orbit-inferred  center 
path  and  the  displacement  anomaly  path  are  virtually  identi¬ 
cal.  Consequently,  it  is  only  necessary  to  display  one  set  of 
calculations.  The  orbit-determined  centers  were  used  here. 

Figure  lOd  illustrates  how  close  the  movements  of  the  two 
observed  and  simulated  rings  are.  There  the  paths  of  the  in¬ 
ferred  centers  along  with  the  displacement  anomaly  path  have 
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been  superimposed.  Except  for  the  beginning  and  end  of  the 
paths,  the  maximum  deviation  is  70  km  with  a  rms  deviation 
of  12  km.  Moreover,  all  paths  show  the  same  general  charac¬ 
teristics.  Both  the  observations  and  simulation  show  a  north¬ 
ward  migration  along  the  isobaths  but  with  some  eastward  or 
looping  component. 


LONGITUDE 
Fig.  lOd 

Fig.  10.  (cont.)  Drifter  paths,  inferred  centers  and  orbit  ellipses  for 
(c)  the  simulated  ring,  and  (d)  center  paths  along  the  continental  slope. 
Arrows  on  the  paths  depict  KMay  intervals,  while  the  ellipses  are 
shown  at  IS-day  intervals. 

The  northward  migration  along  the  slope  is  puzzling,  as 
conventional  theory  [Smith  and  O'Brien,  1983]  would  have 
topographic  beta  drive  the  rings  to  the  south.  Recently,  Smith 
[1986]  has  suggested  that  the  northward  migration  may  be 
due  to  boundary  effects.  An  alternate  explanation  was  offered 
by  Sakamoto  [1986],  who  found  northward  propagating  soli- 
ton  solutions  in  a  two-layer,  nonlinear,  quasi-geostrophic 
model. 

The  swirl  and  translation  velocity  iim?  series  are  shown  in 
Figure  11.  In  the  1980-1981  ring  the  dominant  period  is  15 
days,  almost  exactly  that  found  in  the  midgulf.  For  the  1982- 
1983  ring  the  dominant  period  is  about  11  days,  slightly  less 
than  that  found  in  the  midgulf.  The  simulated  ring  also  shows 


a  decrease  in  the  dominant  period  to  about  15  days.  The 
magnitudes  of  the  observed  ring’s  swirl  velocities  are  about  40 
to  SO  cm/s,  which  is  close  to  that  found  in  the  midgulf.  The 
swirl  velocity  for  the  simulated  ring  is  about  80  cm/s,  substan¬ 
tially  more  than  what  was  obtained  for  the  midgulf.  All  trans¬ 
lation  velocities  indicate  a  fairly  steady  northward  movement 
of  about  4  cm/s. 
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Fig.  11c 

Fig.  11.  Swirl  (left  scale)  and  translation  velocities  (right  scale)  for 

the  (a)  1980-1981,  (6)  1982-1983,  and  (c)  simulated  rings  along  the 
continental  slope. 
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5.  Conclusions 

The  paths,  translation  and  swirl  velocities,  and  orbit  ellipses 
have  been  calculated  from  Lagrangian  observations  for  two 
large  anticyclonic.  Gulf  of  Mexico  rings  and  one  simulated 
nng.  Within  the  capabilities  of  the  data  base  and  the  analysis 
routine,  the  following  are  now  established. 

1.  The  paths  of  the  three  rings  across  the  central  Gulf  of 
Mexico  to  the  continental  slope  off  of  Mexico  are  virtually 
identical.  However,  not  all  Loop  Current  rings  will  follow  this 
path.  See  Lewis  and  Kirwan  [1987]. 

2.  The  intraring  variability  of  translation  and  swirl  veloci¬ 
ties  and  the  inferred  center  paths  for  the  1980-1981  ring  are 
essentially  the  same  as  the  inter-rmg  variability  between  the 
1980-1981  and  1982-1983  rings. 

3.  For  the  observed  rings,  there  was  evidence  of  strong 
ellipticity  developing  as  these  rings  approached  the  conti¬ 
nental  slope.  The  axis  of  orientation  was  approximately 
NE-SW.  No  independent  data  on  ellipticity  were  available  for 
the  1982-1983  ring,  so  it  cannot  be  ruled  out  that  this  ring 
was  so  orientated.  Similar  but  less  pronounced  ellipticity  was 
seen  in  the  simulated  ring. 

4.  Both  observed  rings  showed  significant  distortion  of  the 
drifter  orbits  west  of  94°W.  This  is  attributed  to  interaction  of 
these  rings  with  the  continental  slope  topography  and/or  slope 
circulation  features.  Such  interaction  was  not  observed  in  the 
simulated  ring. 

5.  After  encountering  the  continental  slope  off  of  Mexico, 
both  observed  and  simulated  rings  retained  their  integrity  for 
at  least  several  months  while  they  migrated  to  the  north.  The 
track  of  the  center  of  both  was  very  close  to  the  1200-m 
isobath.  The  observed  swirl  speeds  were  slightly  changed  from 
the  values  obtained  from  the  midgulf.  The  simulated  ring  ex¬ 
hibited  a  substantial  increase  in  swirl  velocity.  For  both  the 
observed  and  simulated  data,  the  translation  speeds  were  com¬ 
parable  to  those  obtained  for  the  midgulf. 

There  are  a  couple  of  broader  implications  of  this  study 
pertinent  to  eddy-resolving  GCM’s.  Both  the  volume  of  data 
as  well  as  the  time  and  space  scale  resolution  of  simulated 
data  available  from  these  models  far  surpasses  the  typical  ob¬ 
servational  data  bases  used  for  verification.  Because  of  the 
different  time  and  space  scales  of  resolution  in  the  simulated 
and  observed  data  bases,  it  is  not  clear  how  to  establish  reli¬ 
able  criteria  for  statistical  comparisons.  This  study  suggests 
that  a  potent  alternative  is  to  compare  model  and  observed 
Lagrangian  kinematics.  Both  simulated  and  observed  La¬ 
grangian  data  sets  provide  comparable  space-time  resolution 
on  the  evolution  of  specific  circulation  features.  Moreover,  the 
assessment  of  the  prediction  of  such  features  is  a  more  strin¬ 
gent  test  of  a  model’s  predictive  capability  than  are  general 
statistics.  The  Lagrangian  data  also  provide  a  means  of  fine- 
tuning  model  parameters  such  as  layer  depth,  density  differ¬ 
ences,  etc.,  so  that  observed  and  simulated  swirl  velocities  can 
be  matched.  As  to  the  HT  model,  this  study  has  documented 
the  interesting  situation  that  the  prediction  improves  in  time, 
at  least  as  far  as  the  movement  of  the  ring  is  concerned. 

The  second  and  related  issue  is  the  use  of  Lagrangian  obser¬ 
vations  for  updating  prognostic  models.  In  character,  such 
data  are  similar  to  sea  surface  topographic  data  derived  from 
satellites  [see  Hurlburt,  1986;  Thompson,  1986;  Kindle,  1986]. 
The  Lagrangian  data  contain  information  on  the  motion  of 
representative  parcels  but  no  information  on  what  is  happen¬ 
ing  nearby.  Satellite  topographic  data  return  information  just 


along  the  satellite  path  but  give  no  information  on  the  actual 
motion  Supplementing  the  latter  data  with  the  former  data  in 
eddy-resolving  general  circulation  models  could  significantly 
upgrade  their  predictive  capabilities. 

Appendix 

The  purpose  here  is  to  provide  some  details  on  the  inver¬ 
sion  of  (3)  and  (4)  to  obtain  the  time  series  for  the  elements  of 
Z^.  The  basis  for  this  inversion  is  a  Taylor  expansion  in  time 
about  the  instant  of  the  velocity  vector.  From  the  left-hand 
side  of  (3)  and  (4)  one  obtains 

u(t)  =  u{tt)  +  M'fttXf  -  t»)  +  -  t»)V2 

+  u'"(.\Xt  -  gV6  (Al) 

lit)  =  iXtJ  -  u'(tjXi  -  ti)  +  v''(t^Xt  -  tk)^/2 

+  ntiXt  -  MV6  (A2) 

for  the  interval  tj  g  t  g  r^*,.  The  primes  represent  time  de¬ 
rivatives.  Now  each  of  the  terms  on  the  right-hand  side  of  (At) 
and  (A2)  can  be  estimated  from  the  velocity  record.  For  exam¬ 
ple  tt(t»)  and  v{t^)  are  merely  the  velocities  at  time  The 
derivatives  can  be  estimated  by  a  variety  of  techniques;  here 
we  have  employed  centered  finite  differences. 

Now  the  analytic  solution,  (5)  and  (6),  can  be  expanded  in  a 
Taylor  series  as  well.  When  this  is  done  and  coefficients  of  the 
appropriate  powers  in  this  expansion  and  (Al)  and  (A2)  are 
equated,  a  system  of  simultaneous  nonlinear  equations  is  ob¬ 
tained  for  each  time  interval.  With  the  subscript  k  suppressed, 
these  are 


X{a  +  d)+  Y{b  -  c)  -I-  2Uj  =  2m 

(A3) 

X{b  +  c)  -f-  Y{d  -  a)  +  IVj.  =  2v 

(A4) 

X[{d  +  a)'  +  b^-  c*]  +  2Yd{b  -  c)  =  4m' 

(A5) 

2Xd{b  +  c)+  y[(d  -  a)^  +  b^  -  c^]  =  4v’ 

(A6) 

Xlid  +  a)’  -1-  (6*  -  c^X3d  +  «)] 

+  Y{b  -  cX}d^  +  a^  +  b^-  c^)  =  8m" 

(A?) 

X{b  +  cX3d^  +  a^  +  b^-  c^) 

+  y[(d  -  a)’  +  (6"  -  c"X3d  -  a)]  =  8m" 

(A8) 

^■[(o^  +  —  c^  +  d*)’  -I-  4d*(<i^  +  b^  —  c^) 

-t-  4ad(o*  +  6*  —  c*  +  d^)] 

+  4Ydlb  -  cXa^  +  b^  -  +  d*)]  =  16u"'  (A9) 

4Xd{b  +  cXfl*  +  b^-c^  +  d*)  -I-  y[(fl^  +  b^-c^  4-  d')^ 

-I- 4d*{a^ -t- 6^ — c’) — 4ad(a^  +  b^—c^+d^)']  =  l 6v"'  (A  10) 

These  eight  equations  are  inverted  at  each  time  step  for  X,  Y, 
Uj,  Yj,  a,  b,  c,  and  d.  To  our  surprise  it  seems  that  this  can  be 
done  analytically.  The  key  to  this  is  the  observation  that  the 
geometric  invariants  of  the  matrix  Af*,  i.e.,  Tr  (M)  and  det  (Af), 
can  be  calculated  from  observations  without  knowing  a,  b,  or 
c,  a  prion.  After  3  years  of  inspection  it  was  seen  that 

det  (Af )  =  2(u"u'"  -  u'"v")l(u‘o”  -  v'uT)  =  Af*  (All) 

Tr  (Af)  =  (u'v'"  —  v'u"')l(u'v"  —  v'u")  =  d  (A12) 

Insertion  of  Af^  and  d,  calculated  from  (All)  and  (A  12), 
respectively,  into  (AS)-{A10)  significantly  simplifies  the  latter 
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expressions.  Considerable  routine  algebra  then  yields 


=  8[u'(M^  -1-  2d^)  -  2u""\/M* 

(A13) 

=  +  2d^)  -  2v""]/M* 

(A  14) 

=  +  K2  +  Ki(d/2)}/N 

(A  15) 

=  + Hi  +  Hi(dl2)}/N 

(A  16) 

=  -  Gj  +  Gi(dl2)]IN 

(A  17) 

where 

N  =  8(u'v"  —  u"v')i 
Gi  =  2(u'^  +  v'^); 

Gj  =  8{u'*  +  v"^); 

Gj  =  16(u'm''  +  v'v''); 

H,  =  2(u'^  - 
Wj  =  8(u"^  -  v"^); 

//j  =  16(u'u"  —  v''v')‘, 

K,  =  2uV; 

Kj  =  8uV': 

A!  3  =  8(u'i)"  4-  u''v'). 

Incidentally,  the  parameter  N  is  directly  related  to  the  angu¬ 
lar  momentum  per  unit  mass  of  the  orbit  L: 

L  =  -4NIM*  (A18) 

The  final  step  in  the  calculations  is  to  determine  the  trans¬ 
lation  velocities  from  (A3)  and  (A4).  This  is  trivial,  since  all 
other  terms  are  now  known. 
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Notes  on  the  Cluster  Method  for  Interpreting  Relative  Motions 

A.  D.  Kirwan,  Jr. 

Department  of  Oceanography,  Old  Dominion  University,  Norfolk,  yirginia 

ll  IS  shown  that  the  so-called  cluster  method  for  determining  the  horizontal  gradients  ol  the  horizontal 
velocity  has  some  basic  limitations  caused  by  restrictive  physical  assumptions.  For  large-area  clusters 
there  is  aliasing  of  these  estimates  by  large-scale  shear,  while  for  small-area  clusters,  the  displacement  of 
the  statistical  cluster  center  from  a  physically  significant  Row  center  can  produce  seriously  aliased  cluster 
model  parameter  estimates.  These  ejects  are  independent  of  typical  experimental  problems  such  as 
measurement  error  and  statistical  reliability.  An  alternative  to  the  cluster  paradigm  is  proposed.  This, 
however,  requires  solving  a  nonlinear  inverse  problem  to  recover  estimates  of  the  horizontal  velocity 
gradients. 


Introduction 

Estimates  of  the  horizontal  spatial  derivatives  of  the  near- 
surface.  horizontal  velocity  are  of  obvious  importance  in  a 
wide  range  of  scales  in  ocean  dynamics.  They  are  notoriously 
difTicult  to  obtain.  Most  estimates  have  come  from  relative 
motion  analysis  of  drifter  clusters.  For  an  older  summary  of 
results,  see  Kirwan  [1975].  Fahrbach  el  al.  [1986]  have  recent¬ 
ly  applied  this  technique  to  drifter  data  from  the  equatorial 
Atlantic. 

Sanderson  el  al  [1988]  have  pointed  out  a  statistical  bias  in 
using  this  technique  if  the  number  of  drifters  in  the  cluster  is 
small  Of  course,  the  statistical  significance  of  model  parame¬ 
ter  estimates,  as  well  as  the  effect  of  measurement  errors  in 
estimating  the  velocity  gradient  tensor,  is  important.  There  is, 
however,  another  problem  with  this  technique.  This  is  the  tacit 
assumption  that  the  center  of  mass  of  the  cluster  corresponds 
to  the  center  of  the  flow  field  about  which  a  Taylor’s  ex¬ 
pansion  of  the  velocity  field  can  be  made. 

The  purpose  here  is  twofold.  The  first  is  to  show  that  with 
the  Cluster  paradigm,  a  deviation  of  the  cluster  center  from  a 
flow  center  can  result  in  serious  biases  in  gradient  estimates. 
The  other  is  to  propose  an  alternative  to  the  cluster  model 
that  avoids  the  center  bias  problem.  The  new  approach  pro¬ 
vides  independent  estimates  of  the  velocity  gradient  estimates; 
however,  it  requires  solving  a  nonlinear  inverse  problem. 

Theory 

The  first  part  of  the  cluster  paradigm  is  given  by 

(r/  -  =  u/  =  £//  -I- 1  +  r^/^  (I) 

Here  A  is  the  time  difference  between  fixes,  r/  is  the  a  compo¬ 
nent  of  the  position  vector  from  the  cluster  center  of  mass  to 
the  drifter  path,  r^/  is  its  original  or  previous  position  vector, 
ff„|,  IS  the  velocity  gradient,  and  UJ’  is  a  steady  mean  flow. 
Figure  la  is  a  cartoon  of  the  general  geometry.  As  has  been 
noted  by  many  previous  investigators,  (1)  can  be  obtained  by  a 
Taylor's  expansion  of  the  velocity  field. 

The  first  term  on  the  right-hand  side  of  (1)  is  a  uniform 
translation  of  the  center,  and  the  second  term  is  a  local  swirl 
induced  by  a  homogeneous  deformation  field.  The  last  term  in 
(1)  IS  a  residual,  the  minimization  of  which  yields  estimates  of 
U/  and  For  the  latter  opieration  it  is  necessary  to  intro- 
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duce  the  cluster  average  of  property  h’’ 

</,<’>  =  Y.  h^/N  (2) 

where  N  is  the  statistical  degrees  of  freedom  in  the  cluster 
average  and  the  sum  is  over  all  drifters  in  the  cluster  [see 
Sanderson  el  al..  1988],  Note  that  by  definition 

<r/>  =  0  (3) 

This  is  the  second  part  of  the  cluster  paradigm.  It  is  stressed 
that  (3)  significantly  constrains  the  basic  assumption  and  va¬ 
lidity  of  a  two-term  Taylor  expansion. 

By  standard  minimization  augmented  with  (3).  the  least 
squares  estimates  of  U/  and  are 

U/  «  <«/>  (4) 

=  I  (5) 

i 

where 

(6) 

is  the  inverse  of  the  statistical  moment  of  inertia  matrix  of  the 
cluster.  Details  and  variations  of  (4)  and  (5)  are  given  by  Molin~ 
ari  and  Kirwan  [1975],  Okubo  and  Ebbesmeyer  [1976],  and 
Okubo  el  al.  [1976], 

This  model  is  critically  dependent  upon  (1)  a  translation 
velocity  Uf  common  to  all  drifters,  (2)  a  velocity  gradient 
tenser  H„i,  common  to  all  drifters,  i.e.,  a  homogeneous  defor¬ 
mation  field,  and  (3)  equation  (3). 

In  order  to  assess  the  role  these  assumptions  play,  consider 
a  more  general  kinematic  flow  model: 

=  K,'  +  Z  -f  r;'VA  (7) 

6 

Here  x/  is  the  position  vector  to  the  drifter  path  from  a 
physically  relevant  local  origin.  This  is  called  the  flow  center 
to  distinguish  it  from  the  cluter  center.  See  Figure  la  for  the 
geometric  relation  between  and  r/.  The  latter  depends 
upon  the  geometric  characteristics  of  the  drifter  array;  the 
former  depends  on  true  flow  field  characteristics.  This  model 
differs  from  (1)  in  that  it  allows  for  a  separate  Taylor's  ex¬ 
pansion  for  each  drifter.  Thus  there  are  distinct  translation 
velocites  y/,  as  well  as  distinct  velocity  gradients  for 
each  drifter.  Also  in  contrast  to  (3) 

<x/>  /  0  (8) 

It  should  be  clear,  however,  that  if  the  three  assumptions  listed 
in  the  preceding  paragraph  are  imposed  on  (7),  it  will  reduce 
to(l). 
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Fig.  1  Cartoon  showing  the  essential  geometry,  (a)  Basic  vectors 
z,  r,  and  x  for  the  present  configuration,  (b)  Displacements  for  case  I. 
This  could  result  from  equal  shear  and  normal  deformation  and  anti- 
cyclonic  rotation.  The  solid  line  indicates  the  present  configuration, 
rnd  the  dashed  line  indicates  the  previous  configuration,  (c)  Displace¬ 
ments  for  case  II  Initially,  ihe  cluster  and  flow  centers  coinade,  but 
in  the  present  (solid  line)  configuration  the  flow  center  has  meved 
away  from  the  cluster  center.  All  panels  depict  the  movement  and 
deformation  of  the  material  lines  connecting  the  drifters  as  well  as  a 
hypothetical  array  initially  symmetne  about  the  flow  center. 


Given  only  cluster  data,  the  observables  are  U/  and  r/. 
We  now  examine  the  effect  of  naively  applying  the  cluster 
paradigm  to  this  more  general  situation.  Taking  the  cluster 
center  as  the  local  origin  and  letting 

z/  =  ’■/  -  W 

be  the  displacement  vector  from  the  cluster  center  to  flow 
center,  it  is  deduced  from  (4),  (5),  (7),  and  (9)  that 

V/  =  <  V/y  +  X  <G^'(r/  -  r/)>  (10) 

» 

=  I  «»,«  K%>  +  I  -  Zs')'-/))  (1 1) 

V  ^ 

These  last  two  equations  show  that  in  general,  the  cluster 
model  parameters  are  weighted  averages  of  the  flow  model.  In 
order  for  (10)  and  (11)  to  reduce  to  (4)  and  (5),  all  three  as¬ 
sumptions  given  above  must  be  made. 

The  following  special  cases  illustrate  the  tenuous  character 
of  'he  cluster  model. 

Case  /.  There  is  no  translation  (V/  sO),  and  the  defor¬ 
mation  field  is  homogeneous;  i.e.,  is  the  same  for  all 
drifters  [G^/  =  G^).  An  example  of  this  situation  would  be 
when  the  drifters  are  embedded  in  a  flow  field  characterized  by 
no  translation  and  homogeneous  vorticity  and  deformation 
fields  (see  Figure  Ih).  It  is  seen  immediately  from  (10)  and  (11) 
that  the  cluster  model  will  yield  bogus  parameter  estimates  of 

=  -  I  G^<z/>  (12) 

b 

(13) 

The  bogus  translation  velocity  has  the  structure  of  a  swirl 


velocity  and  is  of  the  order  of  the  true  velocity  gradient  times 
the  distance  between  the  two  centers.  For  small  area  clusters, 
this  displacement  could  well  be  of  the  cluster  horizontal  scale, 
and  so  the  cluster  determined  translation  velocity  could  be  of 
the  order  of  the  swirl  velocity.  In  the  case  of  the  velocity 
gradient  estimate,  there  is  no  a-priori  reason  to  expect 
<Zz%'’>  to  small;  hence  the  “correction"  term  in  (13)  could 
be  of  the  order  of  the  correct  value.  As  an  aside,  (13)  is  recog¬ 
nized  as  an  extension  of  Pappas’  theorem  from  mechanics. 

Case  11.  The  translation  velocity  varies  from  drifter  to 
drifter;  however,  assumptions  2  and  3  hold  (see  Figure  tc). 
This  case  would  include  the  situation  where  the  drifters  are 
deployed  in  a  jet  such  as  the  Gulf  Stream.  Perhaps,  fortu¬ 
itously,  the  cluster  center  initially  coincides  with  the  flow 
center,  say,  the  axis  of  the  jet.  The  cluster  parameter  estimates 
now  are 

t//  =  <F/>  (U) 

=  +  (15) 

Y 

If  the  flow  center  is  indeed  the  jet  axis,  then  the  cluster 
model  will  yield  a  translation  velocity  for  the  center  somewhat 
less  than  the  maximum  velocity  of  the  jet  as  indicated  by  (14). 
Moreover,  if  the  jet  velocity  is  symmetric  about  the  flow  axis, 
then  the  cluster  gradient  estimates  of  the  cluster  scale  velocity 
gradients  will  include  an  alias  from  the  larger-scale  shear 
across  the  jet  as  exemplified  in  the  inhomogeneous  nature  of 
y/  (the  second  right-hand  side  term  in  (15)).  In  general,  ihe 
cluster  model  may  not  be  able  to  distinguish  between  small- 
scale  and  large-scale  shears. 

These  two  cases  highlight  a  basic  difficulty  of  the  cluster 
method.  If  the  area  covered  by  the  cluster  is  small,  there  is  a 
high  likelihood  that  there  will  be  a  significant  displacement 
between  the  cluster  and  flow  centers.  This  is  where  case  I 
applies.  This  can  be  alleviated  by  increasing  the  area  extent  of 
the  cluster  (with  or  without  additional  drifters),  but  then  the 
large-scale  shear  may  defeat  the  cluster  model,  as  is  indicated 
by  the  analysis  in  case  II. 

Case  HI.  Only  assumption  2  holds.  In  this  case,  the  clus¬ 
ter  parameter  estimates  will  superpose  (12)  and  (14)  for  U/ 
and  (13)  and  (15)  for 

It  is  emphasized  that  the  problems  indicated  above  cannot 
be  rectified  by  increasing  the  number  of  drifters  and/or  the 
cluster  area  and/or  the  position  fix  accuracy  and  frequency. 
The  basic  issue  here  is  the  weakness  of  the  physics  implied  by 
the  approach,  not  the  accuracy  or  statistical  significance  of  the 
measurements. 

Discussion 

The  preceding  analysis  suggests  that  the  cluster  method  for 
estimating  velocity  gradients  should  be  used  with  extreme  cau¬ 
tion.  Other  observations  should  be  used  to  establish  its  appli¬ 
cability  in  specific  experiments. 

There  is  an  alternative  to  the  cluster  paradigm  for  extract¬ 
ing  velocity  gradient  information  from  path  data.  The  alter¬ 
native  makes  use  of  a  kinematic  model  first  proposed  by 
Okubo  [1970]  and  recently  discussed  by  Perry  and  Chong 
[1987].  This  approach  makes  use  of  the  fact  that  (7)  (and 
hence  (1))  can  be  readily  integrated  to  obtain  particle  paths.  If 
the  deformation  field  is  homogeneous  and  stationary,  then  a 
single  particle  path  contains  the  same  information  on  the  de¬ 
formation  as  does  a  hypothetical  cluster  of  nearby  particles. 
The  only  characteristic  difference  in  the  paths  is  their  starting 
positions. 

To  illustrate  how  the  deformation  information  is  contained 
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m  the  path  data,  it  is  necessary  to  calculate  the  path.  Neglect¬ 
ing  the  last  term  in  (7)  and  assuming  that  y/  is  steady  in  time, 
standard  techniques  yield 

w,  "(t)  -  w'.ng  =  f'/ic  -  g  +  x/io  (i6) 
-  g  +  xAO  (17) 

Here  w„'’(t)  is  the  absolute  position  at  time  t,  and  IK’”(to)  “ 
the  initial  position  of  the  flow  origin  at  time  t,  of  drifter  p. 

The  other  quantities  in  (16)  and  (17)  are 

.x/(f)  =  {exp  \_mAt  -  tMXA’"  +  G”)  t 

-  {exp  [m  .  '(t  -  t„)]C,lf ,  '(G'’  -n  +  X^oG , 

(18) 

xA‘)  =  {exp  ImAit  -  giCJC/  Gj/  +  XAI"  -  G')]} 

+  {exp  [m  .  ”(1  - 1.)][  -  Jf .  <-G  j .  <-  +  XAl"  +  G‘^]}/2/' 

(19) 


=  [rr(G,/)±/<-]/2  (20) 

G'’  =  G,/-Gjj'’  (21) 

I"  =  {(G")^  +  (G,/  +  Gj.'O^  -  (G,/  -  (22) 

X/  =  XA‘o)  (23) 


Given  the  path  data,  i.e.,  the  left-hand  side  of  (16)  and  (17) 
for  a  few  intervals  of  time,  the  velocity  gradient  G,/  and  trans¬ 
lation  velocity  V/  can  be  estimated.  This,  however,  is  a  prob¬ 
lem  in  nonlinear  parameter  estimation  and  hence  of  more 
fundamental  mathematical  interest  than  purely  linear  statis¬ 
tical  models.  Kirwan  et  al.  [1988]  have  provided  one  method 
of  obtaining  these  estimates.  Presumably,  more  efficient  meth¬ 
ods  are  available. 

It  is  stressed  that  the  estimates  of  the  velocity  gradient  and 
translation  velocity  apply  only  to  drifter  p.  Other  drifters,  if 
present,  would  supply  independent  estimates.  This  would 
allow  statistical  assessments  of  the  uniformity  of  the  flow  fleld 
parameters. 

A  cautionary  note  on  application  of  the  proposed  technique 
should  be  made.  First,  gradient  estimates  from  this  approach 


are  truly  Lagrangian,  and  so  there  may  be  some  difficulty 
relating  them  to  Eulerian  estimates.  Second,  there  are  certain 
pathological  situations  where  the  deformation  fleld  will  pro¬ 
duce  variability  in  only  one  o[  the  velocity  components.  In 
these  cases  it  may  not  be  possible  to  recover  all  the  gradient 
information.  See  Okubo  [1970]  and  Perry  and  Chong  [1987] 
for  a  complete  discussion  of  pathological  situations  in  critical 
points  in  flow  fields. 
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